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ABSTRACT 


The  NAVSTAR  Global  Positioning  System  (GPS)  has  provided  a  quantum  leap  in 
real-time  autonomous  navigation  capabilities.  NASA’s  Space  Shuttle  will  be  receiving  an 
integrated  GPS  capability  in  the  near  fiiture,  and  the  orbiter  Endeavour  has  been  equipped 
with  a  stand-alone  GPS  receiver.  Although  much  data  is  available  regarding  spacecraft 
GPS  receiver  performance  at  higher  altitudes,  little  information  is  available  for  spacecraft 
at  Shuttle  altitudes  of  approximately  400  km  where  drag  and  gravity  effects  are  more 
pronounced. 

GPS  receiver  navigation  solution  data  from  Shuttle  mission  STS-69  was  made 
available  by  NASA  and  provided  an  opportunity  for  evaluating  GPS  performance  in  low 
Earth  orbit.  NASA  ground  tracking  network  and  Tracking  and  Data  Relay  Satellite 
(TDRS)  data  for  this  mission  provided  a  reference  for  comparison.  Analysis  of  the  data 
was  accomplished  using  Satellite  Tool  Kit  (STK)  for  visualization  and  Matlab  routines  for 
data  comparison. 

GPS  navigation  solutions  were  available  for  approximately  65  percent  of  the  STS- 
69  mission,  and  they  generally  coincided  with  the  reference  track.  Differences  between  the 
GPS  navigation  solution  state  vectors  obtained  using  the  Standard  Positioning  Service 
(SPS)  and  the  referaice  state  vectors  produced  RMS  position  differences  between  the 
data  sets  of  about  1500  m.  One  sigma  position  accuracies  of  54  m  in  the  vertical  direction 
and  approximately  1400  m  in  the  downtrack  direction  were  experienced.  Velocity  vector 
magnitude  differences  during  this  period  were  generally  ±  Im/s,  with  a  RMS  velocity 
difference  of  less  than  9  m/s.  One  sigma  velocity  accuracies  of  approximately  4.2  m/s  in 
the  vertical  direction,  2.3  m/s  in  the  downtrack  direction  and  1.5  m/s  in  the  crosstrack 
direction  were  experienced.  A  firm  conclusion  regarding  Shuttle  GPS  accuracies  could 
not  be  drawn  because  all  sources  of  error  were  not  identified.  Based  on  these  results  GPS 
appears  to  be  an  excellent  navigation  source  for  Shuttle  state  vector  information;  however, 
another  navigation  source  such  as  INS  must  be  present  to  provide  a  check  against 
spurious  data  points  and  periods  of  outage. 
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INTRODUCTION 


The  NAVSTAR  Global  Positioning  System  (GPS)  has  provided  a  quantum  leap  in 
real-time  autonomous  navigation  capabilities.  GPS  technology  has  been  applied  in 
numerous  fields  and  has  now  become  a  method  for  autonomous  spacecraft  navigation. 
Flight  experience  with  the  Topex/Poseidon  spacecraft  has  indicated  that  very  precise 
spacecraft  orbit  information  can  be  obtained  based  on  GPS  measurements.  NASA’s 
Space  Shuttle  will  be  receiving  an  integrated  GPS  capability  in  the  near  future,  and  the 
shuttle  Endeavour  has  been  equipped  with  an  operational  stand-alone  system.  Although 
much  data  is  available  regarding  spacecraft  GPS  receiver  performance  at  higher  altitudes, 
little  information  is  available  for  spacecraft  at  Shuttle  altitudes  of  approximately  400  km 
where  drag  and  gravity  effects  are  more  pronounced. 

GPS  receiver  navigation  solution  data  fi'om  Shuttle  mission  STS-69  has  been  made 
available  by  NASA  and  provides  an  opportunity  for  evaluating  GPS  performance  in  low 
earth  orbit.  NASA  ground  tracking  network  and  Tracking  and  Data  Relay  Satellite 
(TDRS)  data  for  this  mission  was  also  available  for  comparison.  The  objective  of  this 
thesis  was  to  compare  Shuttle  GPS  receiver  navigation  solutions  with  the  NASA  reference 
track.  Data  comparison  routines  written  in  Matlab  were  used  to  analyze  both  sets  of  state 
vectors.  The  visualization  features  of  the  orbit  analysis  software  package  Satellite  Tool 
Kit  (STK)  aided  the  analysis  of  the  data  sets. 

A.  GPS  BASICS 


1.  System  Operation 

GPS  was  developed  by  the  Department  of  Defense  (DoD)  as  a  worldwide  satellite- 
based  radionavigation  system  which  provides  users  with  position,  velocity  and  time 
information.  Using  the  concept  of  ranging,  the  system  measures  the  distance  to  several 
satellites  to  determine  position.  Ranging  is  dependent  on  time  synchronicity  among  the 
receiver  clock  and  the  satellites.  The  satellites  transmit  radio  pulses  in  the  form  of  distinct 
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pseudo-random  sequences  at  specific  known  times.  The  receiving  equipment  recognizes 
these  sequences  and  measures  the  precise  time  at  which  the  pulses  arrive.  Based  on  the 
transmission  and  reception  times,  the  time  it  took  the  pulses  to  travel  fi'om  the  satellite  to 
the  receiver  can  be  determined.  Since  the  speed  of  light  is  constant  and  the  travel  time  of 
the  pulses  is  now  known,  the  range  to  the  satellite  can  be  calculated. 

Determining  the  range  to  a  particular  satellite  places  the  user  on  a  sphere  with  a 
radius  equal  to  this  range.  The  intersection  of  two  spheres  of  position  is  a  circle,  and  the 
intersection  of  three  spheres  is  two  points.  Position  can  be  determined  using  three  satellite 
ranges  because  the  ambiguity  of  the  second  point  can  usually  be  resolved  as  an 
unreasonable  solution.  The  intersection  of  four  spheres  of  position  produces  a  single 
point.  (Herring,  1996,  p.  46) 

Ranging  requires  that  the  receiver’s  clock  be  synchronized  with  the  satellite’s 
clock;  however,  in  order  to  support  receivers  with  less  than  perfect  chronometers,  it  is 
necessary  to  apply  correction  and  this  requires  a  fourth  satellite  as  shown  in  Figure  1.1. 
First,  ranges  to  the  four  satellites  are  calculated  using  the  original  receiver  clock  time  and 
are  called  pseudo-ranges.  The  resultant  four  spheres  of  position  would  intersect  in  a 
single  point  if  there  were  no  clock  error.  The  receiver  clock  error  is  then  determined 
algebraically.  Because  there  is  only  one  value  of  clock  error  for  which  the  four  spheres 
can  be  made  to  intersect  at  a  single  point,  the  user’s  position  solution  is  thus  provided  as 
shown  in  Figure  1.2.  (Herring,  1996,  pp.  46-47) 

Arranged  in  six  orbital  planes  with  four  satellites  in  each  plane,  the  constellation 
of  GPS  satellites  consists  of  twenty-four  satellites.  This  ensures  that  a  minimum  of  four 
satellites  are  always  observable  by  a  user  anywhere  in  the  world.  The  satellites  are  in 
circular  orbits  inclined  at  fifty-five  degrees  with  altitudes  of  20,200  km  and  complete  an 
orbit  every  twelve  hours. 
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UNCORRECTED  SPHERE 
OF  POSITION 

CORRECTED  SPHERE 


Figure  1.1.  Clock  error  in  the  receiver  causes  the  range  to  the  GPS  satellite  to  be 
incorrect.  The  resulting  spheres  of  position  (thick  lines)  do  not  intersect  at  a  single  point. 
By  adjusting  the  receiver  clock  the  ranges  are  corrected  and  meet  at  one  point  (thin  lines). 
In  this  figure  this  appears  to  require  three  satellites;  however,  in  three  dimensions  four 
satellites  are  required.  From  Herring,  1996,  p.  46. 
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Figure  1.2.  Navigation  Using  GPS.  From  NAVSTAR-GPS  Joint  Program  Office,  1987, 
p.  1-17. 
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2.  GPS  Accuracy 


GPS  provides  two  levels  of  service,  the  Standard  Positioning  Service  (SPS)  and 
the  Precise  Positioning  Service  (PPS).  Utilizing  spread  spectrum  techniques,  GPS 
satellites  transmit  on  two  L-band  frequencies,  LI  at  1575.42  MHz  and  L2  at  1227.6  MHz. 
Two  spreading  functions  known  as  C/A-code,  or  Coarse/Acquisition  code,  and  P-code,  or 
Precise  code,  are  employed.  Currently  only  the  C/A-code  signal  is  available  to  all  GPS 
users  and  provides  the  SPS,  while  the  P-code  is  available  to  DoD-designated  users  only 
and  is  necessary  for  the  PPS.  Although  user  position  can  be  determined  using  both  codes, 
only  the  P-code  is  transmitted  on  both  frequencies.  This  enables  users  with  P-code  access 
to  perform  a  dual  frequency  comparison  to  compensate  for  ionospheric  effects  while  C/A 
code  users  must  rely  on  a  less  accurate  model  of  the  ionosphere.  Satellite  ephemeris 
(orbital  element)  data,  almanac  data  for  finding  other  satellites  in  the  constellation, 
atmospheric  propagation  correction  data,  satellite  clock-bias  information,  system  time  and 
satellite  status  information  are  provided  by  the  navigation  message  which  is  transmitted  on 
both  frequencies.  (U.S.  Naval  Observatory,  1996,  pp.  1-2) 

Full  accuracy  is  denied  to  non-DoD  designated  users  through  the  use  of  Selective 
Availability  (SA),  which  modifies  the  navigation  message.  SA  reduces  accuracy  by 
altering  the  satellite  ephemeris  data  and  clock  frequency,  a  procedure  which  is  known  as 
dithering.  Anti-spoofing  protects  agmnst  false  satellite  transmissions  by  encrypting  the  P- 
code,  forming  the  Y-code.  (U.S.  Naval  Observatory,  1996,  p.  2) 

The  SPS  has  a  95  percent  probability  of  providing  horizontd  positioning  accuracy 
within  100  meters,  vertical  positioning  accuracy  within  156  meters  and  timing  accuracy 
within  340  nanoseconds.  (U.  S.  Naval  Observatory,  1996,  p.  1)  Typically,  one  sigma 
positioning  accuracy  is  on  the  order  of  50  meters  horizontally  and  75  meters  vertically 
with  a  velocity  accuracy  of  50  cm/s.  The  specification  for  PPS  requires  positioning 
accuracy  of  16  meters  Spherical  Error  Probable  (SEP)  in  contrast  to  the  SPS  specification 
of  76  meters  SEP  and  a  velocity  accuracy  of  10  cm/s.  A  positioning  accuracy  of  10 
meters  SEP  is  usually  experienced.  (Clynch,  1996)  Sources  of  GPS  inaccuracy  and  their 
impact  are  shown  in  Table  1.1. 
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Component 

Size 

Receiver  Clock 

1000  km 

Selective  Availability 

30  m 

Ionosphere 

1  -  30  m 

Atmosphere 

1  -  6  m 

Orbit  Error 

1-3  m 

Satellite  Clock 

1  -  3  m 

Multipath  -  C/A  Code 

0.5  - 150  m 

Multipath  -  P  Code 

0.1-15  m 

Receiver  Noise 

0.1  -  1  m 

Table  1.1.  GPS  Sources  of  Inaccuracy.  From  Clynch,  1996. 

GPS  system  time  is  maintained  by  the  Composite  Clock  which  incorporates  all 
operational  Monitor  Station  and  satellite  frequency  standards.  This  system  time  is 
referenced  to  the  U.S.  Naval  Observatory’s  Master  Clock  within  a  standard  tolerance  of 
one  microsecond.  Over  the  last  several  years  GPS  system  time  has  maintained  a  tolerance 
within  a  few  hundred  nanoseconds  of  the  Master  Clock.  (U.S.  Naval  Observatory,  1996, 
p.  3) 

B.  SHUTTLE  GPS 

An  integrated  GPS  navigation  capability  has  been  developed  for  the  Shuttle.  The 
scheduled  replacement  of  TACAN  air  navigation  stations  has  been  the  primary  motivating 
factor  in  developing  this  capability.  Currently  TACAN  provides  the  primary  Shuttle  entry 
navigation  aid  from  post  blackout  to  final  approach.  During  final  approach  through 
landing  the  Microwave  Scanning  Beam  Landing  System  (MSBLS)  is  the  primary 
navigation  aid.  Presently,  only  the  orbiter  Endeavour  is  equipped  with  GPS  equipment. 
The  other  orbiters  will  be  equipped  as  they  are  upgraded. 

Two  options  were  considered  for  integrating  GPS  into  the  Shuttle  navigation 
system.  The  first  option  was  to  augment  Shuttle  entry  navigation  filters  to  process  GPS 
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measurements.  The  second  option  was  for  the  Shuttle  navigation  system  to  process  state 
vector  outputs  from  the  GPS  receiver.  Either  option  can  be  extended  to  other  Shuttle 
mission  phases  such  as  ascent  and  on-orbit  operations.  The  second  option  of  processing 
GPS  state  vector  outputs  was  selected  because  it  required  fewer  software  changes  and  had 
the  minimum  impact  on  current  system  operation.  Additionally,  it  fecilitated  the  design 
objective  of  a  phased-development  approach  which  would  permit  evaluation  of  GPS 
receiver  performance  during  flight  before  it  was  used  as  the  primary  entry  navigation  aid. 
(Kachmar  et.  al.,  1993,  pp.  313-317)  With  respect  to  developmental  testing  for  STS-69, 
the  GPS  receiver  navigation  solution  was  not  available  to  the  orbiter  in  real  time  even  as  a 
back  up.  Beginning  with  STS-79,  the  navigation  data  will  be  processed  by  the  orbiter 
computers  and  real-time  state  vectors  will  be  available  during  launch  and  landing. 

Using  the  second  approach,  a  GPS  receiver  state  vector  can  be  used  to  reset  the 
baseline  navigation  state.  This  navigation  state  reset  procedure  is  similar  to  a  ground  state 
vector  uplink  vwth  the  important  difference  that  the  GPS  state  vector  quality  assessment  is 
done  onboard  the  Shuttle.  In  order  to  be  selected  for  guidance,  the  GPS  state  vector  must 
meet  a  predefined  transient  tolerance  level.  Transients  can  be  induced  by  the  Shuttle 
mission  phase  environment,  the  GPS  receiver  switching  between  satellites,  and  periods  of 
satellite  reacquisition.  The  GPS  receiver  can  be  provided  aiding  data  by  the  Shuttle’s 
Inertial  Measurement  Units  (IMU)  to  assist  the  receiver  in  determining  its  position  for 
acquisition  of  the  GPS  satellites.  (Kachmar  et.  al.,  1993,  p.  3 17) 

The  Shuttle  was  designated  as  a  DoD-authorized  user  of  PPS  due  to  the  need  to 
guarantee  a  minimum  level  of  entry  navigation  performance  during  times  of  national 
emergency.  It  was  also  desired  that  current  off-the-shelf  receiver  designs  be  used  that 
complied  with  the  GPS  Joint  Program  Office  Requirements  Document.  (Kachmar  et.  al., 
1993,  p.  317)  The  Collins  3M  receiver  was  selected  for  installation  onboard  Endeavour. 
The  five  chaimel  3M  receiver  tracks  four  satellites  simultaneously  while  reading  the 
navigation  message  for  the  next  satellite  to  be  used  by  the  receiver.  It  is  a  P-code 
receiver  with  position  accuracy  of  15  meters  SEP,  velocity  accuracy  of  10  cm/s,  and 
timing  accuracy  of  40  nanoseconds.  (Collins,  1994)  (Nuss,  e-mail,  1996)  Two  antennas 
are  employed  to  provide  maximum  coverage  for  the  receiver.  The  upper-hemi  and  lower- 
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hemi  antennas  are  located  on  the  top  and  bottom  of  the  orbiter  respectively  as  shown  in 


Figure  1.3. 
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Figure  1.3.  Shuttle  GPS  Antenna  Placement.  From  Carpenter  and  Hinkel,  1995,  p.  6. 


C.  STS-69 

Shuttle  mission  STS-69  lifted  off  from  Pad  39-A  at  the  Kennedy  Space  Center 
(KSC)  on  September  7,  1995,  at  11:09:00.052  a.m.  EDT.  Endeavour’s  crew  of  five 
included:  Commander  David  M.  Walker,  Pilot  Kenneth  D.  Cockrell,  Payload  Commander 
James  S.  Voss,  Mission  Specialist  James  H.  Newman,  and  Mission  Specialist  Michael  L. 
Gemhardt.  The  mission  featured  the  second  flight  of  the  Wake  Shield  Facility  (WSF),  a 
four-meter-diameter  saucer-shaped  free-flying  satellite  designed  to  grow  semiconductor 
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films  in  the  ultra-vacuum  created  in  the  wake  of  the  satellite  as  it  moves  through  space. 
The  Spartan  201  ffee-flying  astronomy  satellite  was  also  deployed  and  retrieved. 
Additionally,  the  crew  tested  assembly  techniques  for  the  international  space  station  and 
tested  thermal  improvements  made  to  space  suits  used  during  space  walks.  On  September 
18,  1995  at  7:37:56  am  EDT,  after  10  days,  20  hours,  28  minutes  and  55  seconds,  STS-69 
landed  at  KSC.  After  a  successful  mission  and  traveling  over  four  and  a  half  million  miles. 
Endeavour  touched  down  on  Runway  33  at  the  KSC  Shuttle  Landing  Facility. 

The  WSF  was  employed  to  conduct  several  GPS  experiments  sponsored  by  the 
Texas  Space  Grant  Consortium  and  was  equipped  with  a  TurboRogue  GPS  receiver 
furnished  by  the  University  Corporation  for  Atmospheric  Research  (UCAR).  University 
of  Texas  at  Austin  researchers  and  Johnson  Space  Center  personnel  conducted  the 
experiments  which  were  aimed  at:  assessing  relative  positioning  using  the  TurboRogue 
receiver  on  the  WSF  and  Endeavour’s  Collins  3M  receiver,  evaluating  precision  orbit 
determination  in  a  low-altitude  environment  (400  km)  where  perturbations  due  to 
atmospheric  drag  and  the  Earth’s  gravity  field  are  more  pronounced  than  for  higher 
altitude  satellites  such  as  TOPEX/POSEIDON,  and  determining  atmospheric  temperature 
profiles  using  GPS  signals  passing  through  the  atmosphere.  (Schutz  et.  al.,  1995,  pp  1-2) 
During  STS-69  Endeavour’s  3M  receiver  operated  in  a  SPS  mode. 
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n.  REFERENCE  FRAMES  AND  SHUTTLE  TRACKING 


A.  REFERENCE  FRAMES 

Several  reference  frames  were  used  when  comparing  GPS  navigation  solutions 
with  the  reference  track.  Shuttle  navigation  and  tracking  sources  utilize  the  Aries  Mean- 
of-1950  Earth-Centered  Inertial  reference  frame,  Greenwich  Mean  Time  and  English  units. 
GPS  navigation  solutions  provide  position  data  in  the  WGS  84  Earth-Centered  Earth- 
Fixed  reference  frame  using  Coordinated  Universal  Time  and  metric  units.  GPS  velocity 
data  is  in  a  Shuttle  Up-East-North  frame  and  in  metric  untis. 

1.  Inertial  Reference  Frames 

The  Aries  Mean-of-1950  (M50)  and  J2000  Earth-Centered  Inertial  (ECI)  reference 
frames  are  not  actually  fixed  with  respect  to  the  mean  position  of  the  stars  in  the  vicinity 
of  the  Sun,  or  inertial  space.  These  inertial  frames  are  defined  by  the  direction  of  the 
Earth’s  axis  of  rotation,  or  celestial  pole,  and  the  direction  from  the  Earth  to  the  Sun  on 
the  first  day  of  spring  which  is  known  as  the  vernal  equinox  or  first  point  of  Aries.  The 
celestial  pole  is  affected  by  gravitational  forces,  primarily  those  exerted  by  the  Sun  and 
Moon.  These  forces  produce  a  quasi-conical  motion  of  the  mean  celestial  pole  around  the 
pole  of  the  ecliptic,  or  axis  of  rotation  of  the  Earth’s  orbit  around  the  Sun.  This  lunisolar 
precessional  motion  has  a  period  of  26,000  years.  The  motion  of  the  actual  celestial  pole 
around  the  mean  pole  is  known  as  nutation  and  has  a  period  of  18.6  years.  (Larson  and 
Wertz,  1992,  pp.  94-95)  (Seidelmann,  1992,  p.  12)  The  International  Astronomical  Union 
(lAU)  has  adopted  conventions  for  calculating  general  precession  (lAU  1976) 
(Seidelmann,  1992,  p.  103)  and  general  nutation  (lAU  1980)  (Seidelmann,  1992,  p.  Ill); 
however,  they  do  not  account  for  transient  pole  wander. 

Motion  of  the  vernal  equinox  subsequently  accompanies  the  motion  of  the  celestial 
pole.  As  a  result,  in  order  to  define  an  ECI  reference  fi'ame  the  date  which  specifies  the 
position  of  the  vernal  equinox  must  be  established.  (Seidelmann,  1992,  p.  12)  The  J2000 
frame  is  a  right-handed  Cartesian  coordinate  system  with  its  origin  at  the  center  of  the 
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Earth  and  its  X  and  Z  axes  pointing  towards  the  mean  vernal  equinox  and  mean  rotation 
axis  of  January  1,  2000.  The  M50  coordinate  system  uses  the  mean  vernal  equinox  and 
mean  rotation  axis  of  1950  as  its  reference  directions.  A  constant  transformation  matrix 
has  been  developed  to  perform  the  conversion  from  M50  to  J2000  coordinates. 
(Seidelmann,  1992,  pp.  184-187) 


2.  WGS  84  Reference  Frame 

The  World  Geodetic  System  1984  (WGS  84)  reference  frame  is  a  Conventional 
Terrestrial  System  (CTS)  whose  origin  is  the  center  of  mass  of  the  Earth.  The  system’s  Z- 
axis  is  parallel  to  the  direction  of  the  Conventional  Terrestrial  Pole  (CTP)  defined  by  the 
Bureau  International  de  I’Heure  (BIH).  The  X-axis  is  the  intersection  of  the  WGS  84 
reference  meridian  plane,  which  is  parallel  to  the  Zero  Meridian  defined  by  the  BIH,  and 
the  CTP’s  equatorial  plane.  The  Y-axis  completes  a  right-handed,  earth-fixed. orthogonal 
coordinate  system.  This  reference  frame  rotates  with  the  Earth  and  assumes  that  the  earth 
is  rotating  at  a  constant  rate  around  a  mean  celestial  pole,  the  CTP.  Variations  from  this 
approximation  cause  the  WGS  84  CTS  to  differ  from  the  true  Instantaneous  Terrestrial 
System  (ITS)  which  is  rotating  around  an  instantaneous  pole,  the  Celestial  Ephemeris  Pole 
(CEP).  The  WGS  84  CTS  can  be  related  to  the  J2000  Conventional  Inertial  System  (CIS) 
by  Equation  2.1. 


CTS  =  [A]  [B]  [C]  [D]  CIS  (2.1) 

The  rotation  matrices  account  for;  [A]  polar  motion,  [B]  sidereal  time,  [C] 
astronomic  nutation,  and  [D]  precession.  Matrices  B,  C  and  D  establish  the  relationship 
between  the  CIS  and  the  ITS.  The  ITS  is  related  to  the  CTS  by  matrix  A  which  provides 
the  relationship  between  the  CEP  and  the  CTP.  (DMA,  1987,  pp.  2-1  -  2-3) 
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3. 


Transformation  to  J2000 


The  Shuttle  navigation  system  and  the  NASA  Shuttle  tracking  sources  produce 
state  vectors  using  the  M50  ECI  reference  frame,  and  GPS  state  vectors  are  provided  in 
the  WGS  84  Earth-Centered  Earth-Fixed  (ECEF)  reference  frame.  As  shown  in  Figure 
2.1,  the  M50  frame  is  fixed  in  inertial  space  while  the  WGS  84  frame  rotates  with  the 
Earth.  In  order  to  compare  state  vectors  from  both  systems  it  is  necessary  to  perform  a 
transformation  between  the  two  reference  frames.  This  transformation  requires  an 
intermediate  step  of  converting  the  data  into  the  J2000  ECI  reference  frame.  For  this 
reason  comparison  of  state  vectors  from  the  two  sources  was  conducted  in  the  J2000 
frame. 


Figure  2.1.  ECI  and  ECEF  Coordinate  Systems.  From  Larson  and  Wertz,  1992,  p.  95. 


4.  Up-East-North  Reference  Frame 

GPS  velocity  data  is  provided  in  the  Up-East-North  (UEN)  reference  frame  which 
has  its  origin  fixed  with  the  spacecraft  as  shown  in  Figure  2.2.  The  UEN  system  can  be 
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visualized  using  the  idea  of  the  Shuttle  flying  with  its  belly  facing  the  Earth  much  like  the 
attitude  an  aircraft  would  assume  while  flying  above  the  surface  of  the  Earth.  The 
Shuttle’s  nose  would  point  along  the  velocity  vector.  Up  would  be  equivalent  to  altitude 
above  the  ground,  East  would  point  straight  ahead  from  the  nose  and  North  would  point 
in  the  direction  of  the  left  wing  tip.  Technically,  the  Up  axis  is  in  the  radial  direction  from 
the  center  of  the  Earth.  The  East  axis  is  in  the  local  orbital  downtrack  direction  and 
parallels  the  velocity  vector.  The  North  axis  is  in  the  crosstrack  direction  completing  a 
right-handed  coordinate  system  and  parallels  the  orbit  normal. 


Z 


RADIAL 


Figure  2.2.  Up-East-North  Reference  Frame. 

In  order  to  perform  the  transformation  from  the  UEN  reference  frame  to  the  WGS 
84  reference  frame  two  Euler  angle  rotations  must  be  performed.  A  rotation  about  the  Y- 
axis  by  +5  followed  by  a  rotation  about  the  Z-axis  by  a  value  of  -y  is  required  as  shown  in 
Equations  2.2, 2.3  and  2.4. 
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This  transformation  is  easily  realizable  in  Fortran  and  was  applied  to  GPS  state 
vector  velocity  components. 


5.  Time 


International  Atomic  Time  (TAI)  compares  several  processes  of  physics  such  as 
cesium  frequency  standards  and  hydrogen  masers  to  form  a  standard  timescale  that  can  be 
used  to  uniquely  identify  the  instance  of  time  at  which  an  event  occurs  on  Earth. 
(Seidelmann,  1992,  p.  2)  Uruversal  Time  (UT)  is  based  upon  the  Earth’s  rotation  with 
respect  to  the  Sun.  Since  the  rotation  of  the  Earth  is  affected  by  irregular  forces,  UT  is 
irregular  with  respect  to  TAI.  The  difference  between  TAI  and  UT  is  increasing  irregularly 
by  approximately  one  second  every  eighteen  months,  and  this  difference  is  always 
maintained  as  an  integral  number  of  seconds  known  as  leap  seconds.  Coordinated 
Universal  Time  (UTC)  is  an  atomic  timescale  that  is  kept  in  close  agreement  with  UT 
through  the  addition  of  these  leap  seconds.  (Seidelmann,  1992,  pp.  6-7)  UTC  is  the 
timescale  used  by  GPS  and  is  provided  in  the  receiver  navigation  solution. 

Greenwich  Mean  Time  (GMT)  is  the  basis  of  civil  time  for  the  United  Kingdom 
and  is  subsequently  related  to  UTC;  however,  in  navigation  terminology  GMT  has  been 
used  as  UT.  As  a  result,  the  term  GMT  is  somewhat  ambiguous.  (Seidelmann,  1992,  p.  7) 
All  Shuttle  tracking  data  uses  GMT  as  its  timescale.  In  this  specific  case,  GMT  is 
equivalent  to  UTC.  (Nikolaidis,  e-mail,  June  1996)  In  order  to  facilitate  the  use  of  the 
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coordinate  frame  transformation  program  and  STK,  all  times  were  converted  to  Mission 
Elapsed  Time  (MET)  in  seconds  which  starts  at  the  time  of  launch. 

B,  SHUTTLE  TRACKING 

NASA  uses  two  primary  sources  for  Shuttle  tracking  data,  ground  sites  using  C- 
band  radar  and  Tracking  and  Data  Relay  Satellites  (TDRS).  Several  of  the  NASA 
tracking  sources  are  listed  in  Table  2.1.  C-band  radar  sites  use  range  and  angle  data  and 
TDRS  uses  doppler  measurements  to  create  Shuttle  state  vectors.  Normally  two  to  four 
state  vectors  are  generated  per  orbit,  although,  more  vectors  are  generated  during  periods 
of  high  activity  such  as  rendezvous,  translational  maneuvers  and  EMU  alignments.  With 
the  exclusion  of  rendezvous,  the  three  sigma  position  accuracy  of  the  state  vectors  is  200 
m  in  the  radial  direction,  450  m  in  the  downtrack  direction  and  200  m  in  the  crosstrack 
direction.  A  three  sigma  velocity  accuracy  of  0.45  m/s  in  the  radial  direction,  0.20  m/s  in 
the  downtrack  direction  and  0.25  m/s  in  the  crosstrack  direction  is  typically  achieved. 
During  quiescent  periods,  three  sigma  accuracies  of  100  m  and  0.30  m/s  radially,  250  m 
and  0.10  m/s  downtrack  and  100m  and  0.15  m/s  crosstrack  may  be  obtained.  During 
rendezvous  and  other  highly  active  periods  errors  on  the  order  of  several  kilometers  may 
be  experienced.  (RSOC,  1996,  pi) 


Source 

Identifier 

East  TDRS 

ESTR 

West  TDRS 

WSTR 

Ascension  Island 

ASCC 

Bermuda 

BDQC 

Kwajalein  Island 

KMTC 

Wallops  Island 

WLRC 

Table  2.1.  NASA  Shuttle  Tracking  Sources. 
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m.  DATA  FILES,  TOOLS  AND  PROCEDURE 


A.  DATA  FILES 


1.  GPS  Navigation  Solutions 

GPS  data  files  from  the  3M  receiver  for  each  flight  day  of  STS-69,  Julian  days  250 
through  260,  were  made  available  by  the  Johnson  Space  Center  (JSC).  There  was  an 
average  of  60  navigation  solution  files  per  day  and  a  total  of  660  files.  Each  file  had  from 
800  to  3300  state  vectors  vwth  fixes  obtained  generally  every  second,  although  every  day 
had  some  period  of  outage  with  no  GPS  data.  Individual  state  vector  entries  contained; 
GPS  system  time;  Coordinated  Universal  Time  (UTC);  position  in  XYZ  WGS  84  ECEF 
coordinates  (m);  velocity  in  UEN  coordinates  (m/s);  position  in  latitude,  longitude  and 
altitude;  pitch,  roll  and  heading  information;  acceleration  in  UEN  coordinates;  and  receiver 
health  and  status  information.  A  t5rpical  GPS  data  file  is  shown  in  Figure  3.1. 

2.  NAVG-11  State  Vectors 

Rockwell  Space  Operations  Company  (RSOC)  provided  a  NAVG-11  Shuttle 
Navigation  Postflight  Product  which  included  ground  navigation  history  vectors.  The 
NAVG-11  product  contains  the  orbiter  solution  vectors  which  are  processed  real-time 
using  C-band  range  and  angle  data  and  TDRS  doppler  data.  Each  entiy  contained  the 
batch  ID,  orbit  number,  batch  number,  UTC  labeled  as  Greenwich  Mean  Time  (GMT), 
and  position  and  velocity  in  XYZ  M50  ECI  coordinates  (ft  and  ft/s).  There  were 
approximately  430  state  vectors  for  the  entire  mission  in  this  file.  Two  to  four  vectors 
were  provided  per  orbit  with  the  bulk  of  the  data  coming  from  TDRS.  (RSOC,  1996,  p. 
1 )  A  portion  of  the  NAVG- 1 1  file  is  shown  in  Figure  3.2. 
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3. 


Propagated  State  Vectors 


RSOC  also  provided  a  second  more  densely  propagated  data  file  with  over  16,000 
state  vectors  for  the  entire  mission.  The  first  column  of  the  propagated  file  is  Mission 
Elapsed  Time  (MET)  in  hours  from  launch  which  was  given  to  be  1995:250:15:09:00.0 
GMT.  The  program  that  produces  the  denser  ephemeris  propagates  fi-om  vector  to  vector 
using  a  Nystrom-Lear  fourth  order  integrator  for  both  position  and  velocity.  Four 
derivative  evaluations  are  required  for  each  60  second  integration  step.  The  program 
takes  into  account  modeled  translational  maneuvers.  Sun  and  Moon  induced  gravitational 
perturbations,  drag  computed  using  a  Jacchia  1970  atmosphere  model,  and  a  7x7 
gravitational  potential  model.  Figure  3.3  shows  a  small  portion  of  the  propagated  data 
file.  (Nikolaidis,  e-mail,  March  1996) 

B.  TOOLS 


1.  File  Transfer  Protocol 

All  analysis  was  performed  on  a  Silicon  Graphics  Indigo  11  workstation.  File 
Transfer  Protocol  (FTP)  provided  the  only  manageable  means  for  retrieving  the  GPS  data 
files  from  the  JSC  server.  The  large  number  and  size  of  the  files  made  the  option  of 
copying  them  to  some  form  of  media  for  transfer  unattractive. 

2.  Unix  Scripts 

Unix  scripts,  sequences  of  Unix  commands  that  can  be  stored  in  files,  were 
employed  in  editing  all  the  data  files  and  were  crucial  in  the  editing  of  the  GPS  data  files. 
Hand  editing  all  the  data  files  was  not  feasible,  and  the  Unix  scripts  proved  to  be  the  only 
way  to  deal  with  the  large  quantity  of  data  files  by  automating  the  procedure  of 
uncompressing  and  compressing  files,  grouping  files  by  day,  editing  files  and  executing 
Fortran  programs.  More  specifically,  the  Stream  Editor  (SED)  command  made  the 


18 


CiPS  UTC  ECEFJC  ECEF_Y  ECEF_Z  VEL_E  VEL  N  VEL  U  LAT  LON  ALT_MSL  ALT_ABS  PITCH  ROLL  HEADING  ACCEL_E 
ACCEL  N  ACCEL_U  CHIFAULT  CHIFREQ  CHIODDE  OHSTATE  CHIHWCH  CHIPRN  CHISNR  CHIJS  CH2FAULT  CH2FREQ 
CH2C0DE  CH2STATE  CH2HWCH  CHBPRN  CH2SNR  CH2JS  CH3FAULT  CH3FREQ  CH3CODE  CH3STATB  CH3HWCH  OHSPRN 
CH3SNR  CH3JS  CH4FAULT  CH4FREQ  CH4CODE  CH4STATE  CH4HWCH  CH4PRN  CH4SNR  CH4JS  CH5FAULT  CH3FREQ 


analysis  of  the  enormous  amount  of  data  possible.  It  provided  an  automated  means  for 
removing  text  headers  and  troublesome  characters  such  as  colons  and  slashes. 


Figure  3.1.  GPS  Navigation  Solution  Data  File. 
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Figure  3.3.  Propagated  Data  File. 


21 


0  181829  -0.17W490E+08  0,6705322E+07  0.9410251E+07  -0, 1  l48541E+0i  -0.2271512E+05  ■0,463576«E404 

0.182107  -0.1794637E+U8  0.6682602E-Kr7  {).9403«)gE+07  •0.1145954E+05  -0.22724WE<K)5  -0.464941 6E404 

0,182384  -0.17957828+08  0,6659872E+07  0.9400952E+07  -<i.ll433fi4E+05  -fl.2273475E+05  -O.4663056E+04 

0.182662  4).1796924B+08  0.6637133E+07  0.939fi282E+07  -0.n40774E+05  -0.2274451E+0.N  -0.4676689E+04 

0.182940  ^>.1798063E+08  0.6614383E+07  0.9391599E-H)7  -0.1 13818 1E-M)5  -fl.2273423E+03  -0.4ri9fl313E-K)4 


3. 


Fortran  Programs 


Fortran  programs  proved  to  be  very  versatile  and  effective  in  dealing  with  the 
varying  data  files.  They  were  used  to  select  the  desired  state  vector  elements  from  the 
input  file  line  entries,  convert  the  input  file  times  to  MET  in  seconds,  and  write  the  new 
data  files.  In  the  case  of  the  GPS  data  files,  Fortran  programs  were  also  used  to  thin  the 
dense  GPS  data  files  by  a  factor  of  ten  and  transform  GPS  UEN  velocity  data  into  WGS 
84  coordinates.  Fortran  routines  also  converted  all  the  NAVG-1 1  and  propagated  data 
from  English  units  to  metric  units. 

4.  Satellite  Tool  Kit 

Satellite  Tool  Kit  (STK)  is  an  orbit  analysis  software  package  developed  by 
Analytical  Graphics  that  employs  an  object  paradigm.  A  Graphical  User  Interface  (GUI) 
is  used  to  manipulate  objects  such  as  scenarios,  vehicles,  facilities,  targets,  and  sensors. 
This  thesis  required  work  with  only  scenarios  and  vehicles.  A  vehicle  is  defined  as  a 
movable  land  sea,  air  or  space  object.  A  scenario  is  a  set  of  objects  that  is  to  be 
visualized.  The  objects  in  the  scenario  are  related  by  time  to  a  map  background  which 
provides  a  geographic  reference  frame.  (Analytical  Graphics,  1996,  pp.  1-2)  STK  Vehicle 
files  for  the  Shuttle  were  created  using  GPS  ephemeris  data  and  the  propagated  file 
ephemeris  data  and  incorporated  into  scenarios  for  visualization. 

5.  ConvertTool 

The  STK  math  utility  convertTool  provided  a  means  for  transforming  data  files 
from  one  reference  frame  to  another.  In  particular  it  permitted  files  to  be  transformed 
from  M50  to  J2000,  J2000  to  M50,  ECEF  to  J2000,  and  J2000  to  ECEF.  The  program 
accounts  for  general  precession  according  to  lAU  1976  and  general  nutation  according  to 
lAU  1980.  It  does  not  account  for  the  Earth’s  irregular  rotation  rate  or  pole  wander. 
(Woodbum,  e-mail,  1996) 
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6. 


Matlab 


Matlab  is  an  interactive  numerical  computation,  data  anal5ns  and  plotting  software 
package.  The  basic  Matlab  data  element  is  a  matrix,  and  it  uses  matrix  manipulation  to 
perform  many  numeric  calculations  in  a  fraction  of  the  time  it  would  take  to  write  and  run 
a  program  in  a  language  such  as  Fortran  which  would  perform  the  same  fiinction.  Matlab 
allows  the  execution  of  sequences  of  commands  that  can  be  stored  in  files  called  M-files. 
Matlab  was  used  extensively  to  compare  GPS  data  and  Shuttle  reference  track  data. 

C.  PROCEDURE 

Upon  undertaking  a  project  with  such  a  large  amount  of  data,  procedures  needed 
to  be  established  to  facilitate  the  transfer  and  manipulation  of  the  vast  data  sets.  From  the 
beginning,  using  File  Transfer  Protocol  to  retrieve  the  data,  to  the  end,  using  Matlab  to 
perform  major  coordinate  transformations,  the  emphasis  was  on  using  computer  programs 
and  tools  to  simplify  all  aspects  of  working  with  the  data.  The  procedure  and  tools 
utilized  in  working  with  the  data  is  outlined  in  the  flow  chart  in  Figure  3.4. 

1.  Obtain  Data  Files 

File  Transfer  Protocol  (FTP)  was  employed  to  retrieve  the  660  GPS  data  files  from 
the  JSC  server.  The  more  manageable  NAVG-1 1  and  propagated  files  were  obtained  on 
floppy  disk  and  8  mm  tape,  respectively.  Grouping  the  files  by  flight  day  provided  the 
most  efficient  way  of  handling  the  data.  The  best  results  were  achieved  by  plotting  each 
day’s  worth  of  data  in  both  STK  and  Matlab. 

2.  Edit  GPS  Files  with  Unix  Scripts  and  Fortran  Program 
a,  Automaiion  with  Unix  Scripts 

The  next  step  was  to  edit  and  manipulate  the  files  using  Unk  scripts  and 
Fortran  programs.  The  following  Unix  script,  “thinfiles,”  was  invoked  to  execute  the 
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scripts  and  the  Fortran  program  that  would  manipulate  all  660  GPS  data  files  for  flight 
days  250  through  260. 
thinfiles 


fileOscript 
file  1  script 
file2script 
files  script 
file4script 
files  script 
file6script 
fileTscript 
fileSscript 
filePscript 
file60script 

Each  file#script  in  “thinfiles”  executed  the  “edit  1  script”  script  for  all  the  GPS  data  files  in 

each  flight  day  where  #  was  the  last  digit  in  the  JuUan  date. 
filel  script 

editlscript  251001 
edit  1  script  251002 
editlscript  251003 
editlscript  251004 


editlscript  251086 

The  “editlscript”  Unix  script  was  the  most  important  script.  It  uncompressed  and 
compressed  the  raw  data  files,  ran  the  stream  editor  command,  executed  the  Fortran 
program,  invoked  convertTooI  and  created  the  new  edited  data  files. 


editlscript 


uncompress  Sl.sol 

sed  -f  sedlscript  Sl.sol  >  dat.edt 

a.  out 

cp  out.ecf  Sl.ecf 
convertTooI  <out.ecf>  $l.j20 
rm  dat.edt 
rm  out.ecf 
compress  Sl.sol 


-  uncompresses  raw  data  files 

-  runs  stream  editor  and  creates  data  file 

-  executes  Fortran  program 

-  copies  Fortran  program  output  to  new  file 

-  invokes  convertTooI  program 

-  removes  intermediary  data  file 

-  removes  intermediary  data  file 

-  compresses  raw  datafile  for  storage 
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The  stream  editor  command  eliminated  countless  hours  of  hand  editing  for  the  raw 
GPS  data  files  by  providing  an  automated  way  of  removing  unwanted  file  header  lines  and 
troublesome  colons  and  slashes.  The  script  “sedl  script”  is  called  by  the  SED  command 
and  coirtains  the  editing  commands  that  were  desired  to  be  performed. 


sedl  script 

/GPS/d  -  deletes  file  header  line  which  contains  the  character  string  GPS 

s/;/  /g  -  replaces  all  colons  with  spaces 

s!/!  !g  -  replaces  all  slashes  with  spaces 

b.  File  Editing  with  Fortran  Program 

The  Fortran  program  “thingps”  removed  the  desired  state  vector  elements 

from  the  edited  input  file,  converted  time  to  MET  in  seconds,  transformed  GPS  UEN 

velocity  data  into  WGS  84  coordinates,  created  the  new  data  files,  wrote  file  header 

information  required  to  execute  convertTool  and  thinned  the  GPS  data  files  by  a  factor  of 

ten.  An  example  of  the  output  from  “thingps”  is  shown  in  Figure  3.5. 

thingps.f 
program  thingps 

implicit  none  ■  declare  variables 

double  precision  x,y,z,velx,vely,velz,velup,veleast,velnorth 
double  precision  gps 
integer  yr,doy,hr,min,sec,mo,day 

integer  year,dayoyear,month,dayomon,hour,minute,second 
integer  tmission,secs 
integer  i,count,il 

open(5,file=='dat.edt',status- old')  -  open  edited  datafile 

rewind(5) 

open(6,file='out.ecf,status='unknown’)  -  open  new  datafile 

rewind(6) 


-  data  required  for  convertTool 
header  and  MET  conversion 


year=1995 
dayoyear=250 
month=09 
dayomon=07 
hour=15 
minute=9 
second=0 

tmission=dayoyear*24*3600+hour*3600+minute*60+second 
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-  comertTool  header  information 


write(*,l)  3 

1  fonnat(il) 
write(*,2)  year,month,dayomon,hour,minute,second 

2  fonnat(i4,5(5x,i2)) 

count=0 

do  100  i=l,90000 

read  (5,*,end=200)  gps,yr,doy,hr,niin,sec,x,y,z,  -  read  desired  data 

%  veleast,velnorth,velup 

secs=(doy*24*3600+hr*3600+inin*60+sec)-tmission  -MET conversion 

calluen(x,y,z,velup,veleast,velnorth,velx,vely,velz) 

il  =  i/10*10  -  thin  data  files  by  factor  of  ten 

if  (i.eq.il)  then 

write(6,3)  secs,x,y,z,velx,vely,velz 

end  if 

count=count+l 
100  continue 

200  continue 

close(5) 
cIose(6) 

3  format(i7,6(2x,fl6.6)) 
end 

subroutine  uen(x,y,z,velup,veleast,velnorth,velx,vely,veIz)  -  velocity  conversion 
implicit  none 

double  precision  x,y,z,velup,veleast,velnorth 
double  precision  R,gamma,d,delta,vebc,vely,velz 

R=(x**2+y**2+z**2)**0.5  -  elements  required for  conversion 

gamma=atan2(y,x)  {see  UEN  reference  frame  section) 

d=(x**2+Y**2)**0.5 
delta=atan2(z,d) 

velx=cos(gamma)*cos(delta)*velup-sin(gamma)*veleast-cos(gamma)  -velocity 
%  *sin(delta)*velnorth  components 

veIy=sin(gamma)*cos(delta)*velup+cos(ganima)*veleast- 
%  sin(gamma)*sin(delta)*velnorth 

velz=sin(delta)*velup+cos(delta)*velnorth 

return 

end 
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a  Reference  Frame  Transformation  Using  ConvertTool 
The  ConvertTool  program  transformed  the  Fortran  output  file  “out.ecf’ 
from  ECEF  coordinates  to  J2000  coordinates  producing  a  file  with  the  extension  “.j20”. 
The  “out.ecf’  header  “3”  corresponds  to  this  conversion  to  J2000,  and  a  “1”  would 
correspond  to  a  conversion  fi'om  M50  coordinates  to  J2000  coordinates.  ConvertTool 
also  requires  a  start  time,  or  epoch,  in  the  file  header  when  converting  to  or  fi-om  ECEF 
coordinates.  This  start  time  was  selected  to  be  the  time  of  launch  in  order  to  comply  with 
the  MET  convention.  Figure  3.6  shows  an  example  output  file  from  convertTool  in  J2000 
coordinates. 


3.  Group  GPS  FUes  by  Day 


After  editing  the  raw  GPS  data  files  and  creating  the  new  files  with  time  in  MET 
and  position  and  velocity  in  J2000  XYZ  coordinates,  the  Unix  script  “daysscript”  was 
utilized  to  group  the  data  files  by  day.  First,  empty  files  were  created  for  each  day  and 
labeled  “gpsday#,”  where  #  was  equal  to  the  last  digit  of  the  Julian  date.  Next, 
“daysscript”  moved  the  empty  files  to  a  temporary  file,  “gps.j20,”  for  editing  by  the 
“append#script.”  The  “Append#scipt”  invoked  the  script  “dataappend”  to  append 
sequentially  all  the  files  for  a  flight  day  to  create  a  single  file.  Finally,  the  edited  files  were 
moved  back  to  the  appropriate  “gpsday#”  file. 


mv  gpsday0.j20  gps.j20 
appendOscript 
mv  gps.j20  gpsday0.j20 
mv  gpsday l.j20  gps.j20 
append  1  script 
mv  gps.j20  gpsday  l.j20 


dataappend  251001 
dataappend  251002 
dataappend  251003 


dataappend  251086 


mv  gpsday60.j20  gps.j20 

append60script 

mv  gps.j20  gpsday60.j20 


cp  $l.j20  data.j20 


cat  data.j20  »  gps.j20 
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210632.000000 -6.5927172845e+06  -8.4299477392eH>5  1.1831728757&+06  2,6972220726c+02  -6.8847722257e+03  -3.4078667026e+03 
210642,000000  -6,58962 17468e+06  -9.1179043947^5  1. 14902 13447e+06  3.5515780257e+02  ■6.8734399550c«)3  -3.4230278559i>+K)3 
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4. 


Edit  Propagated  Data  File 


The  propagated  data  file  required  much  less  manipulation.  A  Fortran  program, 
“edtprop,”  was  used  to  add  the  header  for  convertTool,  convert  the  input  file  times  to 
MET  in  seconds  and  convert  the  data  from  English  units  to  metric  units.  A  copy  of  the 
program  is  in  Appendix  A.  ConvertTool  was  then  used  to  transform  the  file  from  M50 
coordinates  to  J2000  coordinates.  Finally,  the  file  of  over  16,000  state  vectors  was  hand 
edited  to  break  it  up  by  flight  day. 

5.  Create  STK  Vehicle  Files 

STK  Vehicle  files  for  each  flight  day  were  created  using  the  Shuttle  STK  vehicle 
library  file  as  a  guide.  A  vehicle  file  was  created  for  each  day  fi'om  each  data  source.  GPS 
ephemeris  data  and  the  propagated  file  ephemeris  data  were  each  inserted  into  a  model 
Shuttle  vehicle  file  in  ECI J2000  coordinates.  A  scenario  was  created  for  every  flight  day 
which  included  a  vehicle  generated  from  GPS  ephemeris  data  and  one  generated  from  the 
propagated  ephemeris  data.  An  example  vehicle  file  is  shown  in  Appendix  C.  The 
scenarios  were  animated  so  that  both  vehicles  could  be  viewed  simultaneously  to  better 
visualize  the  similarities  and  differences  between  the  tracks. 

6.  Matlab  Comparison 

Matlab  was  employed  for  a  detailed  comparison  of  the  GPS  and  propagated  data 
files  for  each  day.  The  Matlab  M-file,  “diffday#.m”  was  executed  to  plot  J2000  XYZ 
position  and  velocity  vector  components  and  the  position  and  velocity  vector  magnitudes 
for  both  data  sets.  A  plot  of  the  data  points  in  three  dimensions  was  also  obtained  using 
“orbplt#.m”  to  assist  in  visualizing  the  orbit  that  was  generated  by  the  GPS  data. 
Appendices  D  and  E  contain  examples  of  these  M-files. 

Since  the  propagated  data  and  GPS  data  could  not  be  compared  point  for  point,  it 
was  difficult  to  obtain  hard  numbers  for  comparison  of  GPS  and  reference  track  data.  In 
an  effort  to  match  the  data  sets  point  for  point,  the  process  of  editing  the  GPS  data  files 
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for  days  251  and  252  was  repeated  with  the  exception  of  thinning  the  data  by  a  factor  of 
ten.  The  NAVG-11  file  was  hand  edited  to  remove  headers  and  unnecessary  columns, 
and  then  run  through  the  stream  editor  to  remove  colons.  The  Fortran  program  “edtsv” 
shown  in  Appendix  B  was  then  executed  to  convert  state  vector  times  to  MET  and 
English  units  to  metric  units.  ConvertTool  was  invoked  to  transform  the  file  from  M50 
coordinates  to  J2000  coordinates.  GPS  and  NAVG-11  state  vectors  were  then  hand 
matched  for  days  251  and  252.  A  Matlab  M-file  was  used  to  compare  position  and 
velocity  vector  components,  position  and  velocity  vector  magmtudes,  and  Root  Mean 
Square  (RMS)  differences  between  the  GPS  data  and  the  NAVG-11  data.  A  copy  of  the 
Matlab  script  “diffinat#.m”  is  shown  is  Appendix  F. 
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IV.  RESULTS 


GPS  data  was  available  for  every  flight  day  of  STS-69;  however,  GPS  data  was 
not  available  for  the  entire  duration  of  each  flight  day.  The  percentage  of  each  day  for 
which  GPS  data  was  available  for  comparison  with  the  propagated  data  varied  greatly 
fi-om  day  to  day.  Day  258  was  covered  most  sparsely  by  GPS  with  data  being  available 
only  12  pecent  of  the  day.  In  contrast,  GPS  data  was  available  for  93  percent  of  day  255. 
The  average  percentage  of  the  day  during  the  mission  for  which  GPS  data  was  available 
was  65  percent.  These  results  are  shown  in  Figure  4. 1 . 

100 
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^  70 
§  60 
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S,  40 
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251  252  253  254  255  256  257  258  259  260 
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Figure  4.1.  Percentage  of  each  day  for  which  GPS  data  was  available. 

A.  MATLAB  COMPARISON  RESULTS 


1.  Day  251 

Day  251  data  was  typical  of  the  results  achieved  during  the  mission.  Figure  4.2  at 
the  end  of  the  section  displays  how  dense  the  GPS  data  is  relative  to  the  propagated  file 
reference  track.  Of  note  is  the  fact  that  the  GPS  data,  displayed  as  asterisks,  has  already 
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been  thinned  by  a  factor  of  ten.  The  propagated  file  is  derived  from  the  original  state 
vectors  which  come  from  several  sources  as  shown  in  Figure  4.3.  The  sources  displayed 
trends  when  compared  with  the  GPS  data,  and  these  trends  appeared  to  be  diurnal  in  some 
cases.  The  list  of  tracking  source  identifiers  is  shown  in  Table  4. 1 . 


Source 

Identifier 

East  TORS 

ESTR 

West  TDRS 

WSTR 

Ascension  Island 

ASCC 

Bermuda 

BDQC 

Kwajalein  Island 

KMTC 

Wallops  Island 

WLRC 

Table  4.1.  NASA  Shuttle  Tracking  Sources. 


A  plot  of  the  position  vector  XYZ  components  in  J2000  coordinates  for  day  25 1 
shown  in  Figure  4.4  indicates  that  the  GPS  data  represented  by  the  heavy  line  follows  very 
closely  with  the  reference  track  represented  by  the  thinner  line.  A  two  hour  gap  in  GPS 
data  is  seen.  This  is  typical  of  the  gaps  in  GPS  data  that  were  observed  every  day.  A  plot 
of  the  velocity  vector  in  XYZ  components  in  J2000  coordinates  in  Figure  4.5  also  shows 
this  gap  in  GPS  data  and  some  spurious  GPS  data  points.  Both  plots  display  sinusoidal 
patterns  for  each  of  the  components  which  are  associated  with  an  inclined  circular  orbit. 

The  position  vector  magnitude  is  equivalent  to  the  spacecraft’s  distance  from  the 
center  of  the  Earth,  or  radius  of  its  orbit,  as  shown  in  figure  4.6.  Plotting  of  the  position 
vector  magnitude  in  Figure  4.7  shows  that  the  radius  is  fairly  constant  as  would  be 
expected  for  the  Shuttle’s  circular  orbit.  This  plot  of  vehicle  altitude  was  very  effective  in 
indicating  when  vehicle  maneuvers  occurred.  The  velocity  vector  magnitude  plot  also 
indicates  that  the  velocity  vector  was  quite  constant  as  was  expected  for  a  circular  orbit. 
Both  the  position  and  velocity  vector  magnitude  plots  show  that  spurious  GPS  data  points 
were  present. 
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The  three  dimensional  plot  of  the  GPS  data  for  day  251  in  J2000  coordinates 
shown  in  Figure  4.8  displays  how  well  the  GPS  data  captured  the  Shuttle’s  circular  orbit. 
The  fact  that  the  orbit  is  inclined,  and  that  the  GPS  data  is  extremely  dense  is  reflected  in 
this  plot.  Erroneous  data  points  are  also  visible. 

Point  for  point  comparison  of  the  GPS  position  vector  magnitudes  with  the 
NAVG-1 1  position  vector  magnitudes  in  Figure  4.9  shows  an  altitude  difference  of  ±150 
m.  Figure  4.10  shows  that  the  velocity  vector  magnitude  difference  between  the  two 
data  sets  was  generally  ±  1  m/s  and  did  not  exceed  4  m/s.  The  RMS  position  difference 
between  the  GPS  data  and  NAVG-1 1  data  displayed  in  Figure  4. 1 1  was  primarily  between 
1400  and  1550  m,  and  the  RMS  velocity  difference  displayed  in  Figure  4.12  did  not 
exceed  9  m/s. 


2.  Typical  Data 

The  position  and  velocity  vector  magnitude  plots  for  day  252  shown  in  Figure 
4.13  display  a  sinusoidal  pattern  associated  with  what  appears  to  be  an  elliptical  orbit. 
Closer  examination  of  the  position  vector  magnitude  plot  scale  shows  that  the  peak  to 
peak  difference  associated  with  apogee  and  perigee  is  approximately  5000  meters.  Gaps 
in  GPS  data  and  spurious  data  points  are  also  experienced.  Day  255  position  and  velocity 
vector  magnitude  plots  displayed  in  Figure  4. 14  are  similar  to  those  for  day  252;  however, 
the  errant  data  points  are  more  severe. 

GPS  data  for  day  258  was  very  sparse  with  only  12  percent  of  the  flight  day 
covered  by  GPS  data.  The  sparse  amount  of  data  allowed  a  detailed  look  at  both  data  sets 
as  shoAvn  in  Figure  4.15  and  displayed  the  imperfect  nature  of  an  actual  vehicle  orbiting 
the  Earth.  Although  there  were  a  few  spurious  data  points,  the  GPS  data  followed  the 
reference  track  so  closely  that  it  was  difficult  to  discern  the  GPS  data  fi'om  the  reference 
track. 
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3. 


Maneuvers 


The  position  vector  magnitude  plot  proved  extremely  useful  in  detecting  Shuttle 
maneuvers.  The  day  253  position  vector  magnitude  plot  shown  in  Figure  4.16  displays  a 
gap  in  GPS  data  just  prior  to  the  maneuver.  A  small  gap  in  GPS  data  is  experienced  at  the 
beginning  of  the  maneuver  which  involves  an  altitude  change  of  10,000  m;  however,  the 
GPS  track  follows  the  reference  track  through  the  rest  of  the  maneuver.  Figure  4.17  also 
displays  an  example  of  GPS  tracking  through  a  maneuver.  This  maneuver  is  more  severe 
involving  a  change  in  altitude  of  60,000  m.  A  loss  of  GPS  data  is  experienced  in  the 
middle  of  the  maneuver  with  GPS  data  being  reacquired  before  completion  of  the 
maneuver. 


4.  Errant  GPS  Data 

There  are  several  examples  of  errant  GPS  data  over  the  course  of  the  eleven  flight 
days  of  STS-69.  One  of  the  best  examples  occurs  on  day  250  and  is  shown  in  Figure  4. 18 
and  Figure  4.19.  Figure  4.18  is  a  three  dimensional  plot  of  the  GPS  data  for  day  250. 
This  plot  shows  a  segment  of  GPS  data  that  fails  to  conform  to  the  Shuttle’s  circular  orbit. 
This  period  of  bad  data  can  be  clearly  identified  in  the  position  vector  component  plot 
shown  in  Figure  4.19.  At  the  start  of  the  day’s  data  the  GPS  data  does  not  follow  the 
reference  track.  Afl;er  a  period  of  approximately  1300  seconds,  the  GPS  data  changes 
drastically  and  matches  the  reference  track. 

GPS  data  for  day  259  has  an  interesting  anomaly  that  is  exhibited  in  the  position 
vector  magnitude  plot  in  Figure  4.20.  A  segment  of  GPS  data  appears  to  be  disjointed 
from  the  rest  of  the  data.  In  the  three  dimensional  plot  of  Figure  4.21  this  period  of  errant 
data  appears  as  a  separate  orbit  at  a  different  inclination.  Several  spurious  data  points  are 
also  present,  particularly  in  the  velocity  vector  magnitude  plot. 

The  position  and  velocity  vector  magnitude  plots  presented  the  most  dramatic 
results  of  errant  data.  The  Day  260  position  and  velocity  vector  magnitude  plot  in  Figure 
4.22  is  an  example  of  this  observation.  The  last  5000  seconds  of  position  and  velocity 
data  for  the  Day  260  GPS  data  varies  significantly  from  the  reference  track.  This  variance 
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Day  251  Mission  Seconds  ^ 

Figure  4.3.  NAVG-11  Data  Sources. 
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Figure  4.4.  Day  251  Position  Vector  Components  in  J2000  Coordinates. 
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Z  Velocity  (m/s) 


Figure  4.5.  Day  251  Velocity  Vector  Components  in  J2000  Coordinates. 
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Figure  4.10.  Day  251  Velocity  Vector  Magnitude  Difference. 


Position  RMS  Difference  (m) 


ESTR 

o 

WSTR 

X 

BDQC 

+ 

WLRC 

• 

KMTC 

46 


Velocity  RMS  Difference  (m/s) 
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Figure  4.19.  Day  250  Position  and  Velocity  Vector  Magnitude. 
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Z  Axis  (m) 


GPS  Orbit  lor  Day  259  In  J2000  Coordinatas 
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Figure  4.21.  Day  259  GPS  Orbit  in  J2000  Coordinates. 
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Day  260  Mission  Seconds 
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Figure  4.22.  Day  260  Position  and  Velocity  Vector  Magnitude 


B.  STK  RESULTS 


STK  proved  to  be  a  valuable  tool  for  visualizing  the  differences  between  the  GPS 
vehicle  tracks  and  the  propagated  vehicle  tracks.  The  GPS  vehicle  is  labeled 
“STS69gpsday#  ”  and  the  propagated  file  vehicle  is  labeled  “STS69propday#”  where  # 
corresponds  to  the  last  digit  of  the  Julian  date.  The  circle  surrounding  the  vehicle 
represents  the  vehicle’s  field  of  view  as  seen  by  a  horizon  sensor  mounted  on  the  Shuttle. 
This  field  of  view  circle  was  helpful  in  discerning  changes  in  altitude  because  it  expanded 
as  the  vehicle’s  altitude  increased. 

The  ground  tracks  displayed  in  Figure  4.23  show  an  example  of  the  extreme 
differences  that  were  occasionally  encountered  between  the  two  data  sets.  This  example 
occurs  on  day  250  and  provides  another  look  at  the  discrepancy  observed  in  the  Matlab 
plots  for  this  day  shown  in  Figure  4.18  and  Figure  4.19.  The  two  tracks  eventually  did 
coincide  as  shown  in  Figure  4.24.  When  the  vehicle  tracks  coincided  in  this  manner  it  was 
difficult  to  distinguish  between  them  because  they  were  plotted  directly  on  top  of  each 
other. 

Figure  4.25  is  more  typical  of  the  STK  plots  obtained.  This  plot  shows  the  entire 
ground  track  for  day  25 1 .  The  GPS  and  propagated  ground  tracks  overlay  directly  on  top 
of  each  other,  and  the  differences  between  the  two  vehicles  are  indistinguishable  because 
they  are  also  plotted  on  top  of  each  other.  Figure  4.26  displays  a  variance  in  the  data  for 
day  251  showing  the  divergent  tracks  of  the  two  vehicles.  Each  day’s  data  analysis 
showed  at  least  one  anomaly  similar  to  this  one.  The  perspective  shot  in  Figure  4.27 
provides  a  better  view  of  the  difference  between  the  two  tracks.  Figure  4.28  shows  a 
similar  error  for  day  253.  The  ejq)anded  field  of  view  for  the  GPS  vehicle  indicates  that 
there  is  a  difference  in  altitude  between  the  two  vehicles. 

A  perspective  view  is  shown  in  Figure  4.29  and  displays  how  well  the  vehicles  and 
orbits  for  day  255  match  each  other.  Figure  4.30  is  a  close-up  of  the  same  picture  and 
exhibits  how  closely  the  vehicles  and  two  sets  of  orbits  overlay  on  top  of  each  other. 
Figure  4.31  further  zooms  in  and  reveals  that  the  vehicles  are  separated,  with  the  GPS 
track  trailing  the  propagated  track  by  a  small  distance. 
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Several  possibilities  exist  as  to  why  the  GPS  data  and  reference  track  data  do  not 
match  directly.  One  possibility  is  that  the  error  is  caused  by  the  transformation  between 
reference  frames  that  does  not  take  into  account  pole  wander  or  the  irregular  rotation  rate 
of  the  Earth.  Another  possibility  is  a  spurious  data  point  from  the  GPS  receiver.  While 
these  factors  may  contribute  to  the  difference  between  the  data  sets,  the  largest  source  of 
error  is  most  probably  the  fact  that  the  files  for  the  GPS  track  and  the  files  for  the 
reference  track  each  had  times  that  were  rounded  to  the  nearest  whole  second.  With  the 
Shuttle  traveling  at  over  7  km/s  this  can  result  in  considerable  position  differences  between 
the  two  data  sets  particularly  in  the  downtrack  direction. 

Figure  4.32  shows  a  GPS  position  data  point  at  the  center  of  a  box  whose  sides 
were  determined  based  on  the  velocity  of  the  Shuttle  in  the  X,  Y  and  Z  directions  for  a 
period  of  one  second  at  a  particular  instant  on  day  251.  This  error  box  is  approximately 
6000  m  long  in  the  X  direction,  2500  m  wide  in  the  Y  direction  and  1200  m  high  in  the  Z 
direction.  As  expected,  the  reference  track  data  point  for  the  same  instant  fell  within  the 
box.  Figure  4.33  shows  a  plot  of  the  differences  between  the  GPS  data  and  the  NAVG-1 1 
state  vectors  for  day  25 1 .  The  majority  of  the  differences  between  the  data  sets  lies  within 
the  box  which  encloses  the  error  which  may  have  been  induced  by  a  full  second  of 
timescale  inaccuracy. 

The  most  extreme  example  of  divergence  between  the  GPS  and  propagated  tracks 
occurred  on  day  256  and  is  shown  in  Figure  4.34.  In  this  plot  the  perspective  view  is  from 
150,000  km  above  the  Earth  and  the  small  disc  at  the  center  is  the  Earth.  Figure  4.35  is 
more  indicative  of  the  differences  between  vehicle  tracks  that  were  encountered.  TypicaUy 
an  error  in  which  the  GPS  track  started  to  lead  or  lag  the  propagated  vehicle  would  occur 
once  per  day. 

The  ground  track  plot  for  day  258  shown  in  Figure  4.36  reflects  the  sparse  data  for 
this  flight  day.  The  perspective  view  in  Figure  4.37  shows  how  closely  the  two  vehicle 
tracks  match.  The  close-up  view  in  Figure  4.38  also  shows  that  the  GPS  track  coincides 
with  the  propagated  track.  Figure  4.39  reveals  that  the  two  tracks  are  right  on  top  of  each 
other.  The  GPS  orbit  in  this  plot  appears  very  smooth  while  the  propagated  orbit  appears 
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to  be  segmented  and  not  as  smooth  as  the  GPS  orbit.  This  is  expected  since  the  GPS  orbit 

is  composed  of  considerably  more  data  points. 

The  plot  for  day  259  shown  in  Figure  4.40  indicates  some  discrepancies  between 
the  GPS  and  propagated  vehicle  ground  tracks.  It  shows  two  independent  vehicles  which 
indicates  that  the  vehicle  tracks  diverge.  Figure  4.41  further  shows  that  the  GPS  data  has 
produced  a  vehicle  in  a  different  orbit  at  a  different  inclination.  This  plot  matches  the 
Matlab  three  dimensional  GPS  orbit  plot  for  this  day  shown  in  Figure  4.21. 
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Figure  4.24.  Day  250  STK  Coincident  Shuttle  Vehicle  Tracks. 
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Figure  4.25.  Day  251  STK  Coincident  Shuttle  Vehicle  Tracks. 
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Figure  4.26.  Day  25 1  STK  Divergent  Shuttle  Vehicle  Tracks. 
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Figure  4.28.  Day  253  STK  Divergent  Shuttle  Vehicle  Tracks. 
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Figure  4.31.  Day  255  STK  Perspective  View  Close-up  of  Shuttle  Vehicle  Tracks. 
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Day  251  GPS  and  Reference  Positions  in  J2000  Coordinates 
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Figure  4.32.  Day  25 1  GPS  and  Reference  Positions  in  J2000  Coordinates 
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Figure  4.34.  Day  256  STK  Perspective  View  of  Divergent  Shuttle  Vehicle  Tracks. 
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Figure  4.38.  Day  258  STK  Perspective  View  Close-up  of  Shuttle  Vehicle  Tracks. 
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Figure  4.40.  Day  259  STK  Divergent  Shuttle  Vehicle  Tracks. 
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c. 


SHUTTLE  UEN  COMPARISON  RESULTS 


1.  UEN  transformation 

In  an  effort  to  better  classify  the  state  vector  errors  encountered  as  errors  in  the 
crosstrack  direction,  in  the  downtrack  direction,  or  in  altitude,  the  GPS  and  NAVG-11 
state  vectors  with  matching  times  were  transformed  from  J2000  coordinates  to  the 
Shuttle’s  UEN  frame.  In  order  to  perform  the  transformation  from  the  J2000  reference 
frame  to  the  Shuttle  UEN  reference  frame  as  shown  in  Figure  4.42,  two  Euler  angle 
rotations  must  be  performed.  A  rotation  about  the  Z-axis  by  a  value  of  y  followed  by  a 
rotation  about  the  Y-axis  by  -6  are  required  as  shown  in  Equations  4.1,  4.2  and  4.3. 


Z 


X 

Figure  4.42.  Up-East-North  Reference  Frame. 
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2.  UEN  results 

This  transformation  was  implemented  using  Matlab  and  is  available  in  Appendix  G. 
The  following  results  shown  in  Tables  4.2  and  4.3  and  Figures  4.43,  4.44,  4.45, 4.49,  4.50 
and  4.51  at  the  end  of  the  chapter  were  obtained  from  the  day  251  and  day  252  state 
vector  position  differences. 


Vertical  (m) 

Downtrack  (m) 

Crosstrack  (m) 

Average 

12.1 

5.5 

0 

Median 

13.2 

85.7 

0 

Standard  Deviation 

53.8 

1461.0 

0 

Table  4.2.  Day  251  State  Vector  Position  differences. 


Vertical  (m) 

Downtrack  (m) 

Crosstrack  (m) 

Average 

2.4 

-32.7 

0 

Median 

8.0 

-140.0 

0 

Standard  Deviation 

53.8 

1348.9 

0 

Table  4.3.  Day  252  State  Vector  Position  differences. 
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The  vertical  position  differences  shown  in  Figure  4.43  and  4.49  were  less  than  the 
expected  one  sigma  value  of  75  m  for  SPS.  The  values  obtained  for  vertical  position 
difference  matched  the  values  obtained  from  the  position  vector  magnitude  (altitude) 
difference  plot  when  the  data  sets  were  analyzed  in  the  J2000  reference  frame.  These 
values  are  within  the  expected  SPS  156  m  vertical  positioning  accuracy. 

Although  the  average  downtrack  position  difference  was  fairly  small  for  both  days, 
the  extremely  large  standard  deviation  indicates  that  significant  downtrack  errors  were 
encountered.  When  traveling  at  7  km/s  a  downtrack  error  of  1400  m  corresponds  to  a  0.2 
s  error  in  time.  It  was  expected  that  a  constant  position  error  in  the  downtrack  direction 
would  be  observed.  The  downtrack  position  differences  displayed  in  Figure  4.44  and 
Figure  4.50  varied  approximately  +1400  m,  and  no  constant  offset  was  discovered.  These 
results  far  exceeded  the  SPS  one  sigma  horizontal  positioning  error  of  50  m. 

Crosstrack  errors  were  all  less  than  1x10'®  m  for  both  days  as  shown  in  Figure  4.45 
and  Figure  4.51.  This  was  beyond  the  precision  of  the  data,  and  makes  the  differences 
essentially  zero.  These  results  were  unexpected  and  did  not  match  the  one  sigma 
horizontal  position  error  of  50  m  for  SPS.  Due  to  unidentified  sources  of  error,  no  firm 
conclusions  can  be  drawn;  however,  the  downtrack  and  crosstrack  errors  seem  to  indicate 
that  most  of  the  position  error  is  in  the  downtrack  direction. 


Vertical  (m/s) 

Downtrack  (m/s) 

Crosstrack  (m/s) 

Average 

0.29 

0.83 

-0.16 

Median 

0.70 

0.11 

-0.13 

Standard  Deviation 

4.18 

2.15 

0.96 

Table  4.4.  Day  251  State  Vector  Velocity  differences. 


Vertical  (m/s) 

Downtrack  (m/s) 

Crosstrack  (m/s) 

Average 

0.58 

1.13 

-0.06 

Median 

0.38 

0.39 

-0.15 

Standard  Deviation 

4.08 

2.34 

1.53 

Table  4.5.  Day  252  State  Vector  Velocity  differences. 


81 


The  velocity  results  from  the  comparison  of  state  vectors  for  day  251  and  day  252 
are  shown  in  Tables  4.4  and  4.5  and  Figures  4.46,  4.47,  4.48,  4.52,  4.53  and  4.54.  Both 
days’  results  are  very  similar.  In  general  the  average  values  for  vertical  and  dovratrack 
velocity  differences  exceeded  the  expected  SPS  one  sigma  velocity  accuracy  of  0.5  m/s, 
and  the  median  values  were  close  to  this  value.  The  crosstrack  velocity  data  is  well  within 
the  expected  SPS  one  sigma  velocity  accuracy. 


Figure  4.43.  Day  25 1  Vertical  Position  Difference. 
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Figure  4.47.  Day  251  Downtrack  Velocity  Difference. 
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Figure  4.49.  Day  252  Vertical  Position  Difference. 
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Figure  4.50.  Day  252  Downtrack  Position  Difference 
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Figure  4.51.  Day  252  Crosstrack  Position  Difference. 
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Vetical  Velocity  Difference  (m/s) 
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Figure  4.53.  Day  252  Downtrack  Velocity  Difference. 


92 


6  1.7  1.8  1.9  2  2.1 

Day  252  Mission  Seconds  .jqS 

Figure  4.54.  Day  252  Crosstrack  Velocity  Difference. 
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V. 


CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

GPS  navigation  solutions  were  available  for  approximately  65  percent  of  the  STS- 
69  mission,  and  they  generally  coincided  with  the  reference  track.  Based  on  vehicle  track 
visualization  using  STK  and  based  on  conversion  from  J2000  coordinates  to  a  spacecraft 
local  UEN  reference  fi^e,  differences  between  the  GPS  data  and  the  reference  data 
appear  to  have  occurred  predominantly  in  the  dowmtrack  direction.  The  fact  that  the  two 
data  sets  were  rounded  to  the  nearest  whole  second  probably  impacted  the  downtrack 
differences  between  them  significantly. 

State  vector  differences  between  the  GPS  navigation  solutions  obtained  using  the 
Standard  Positioning  Service  and  the  reference  NAVG-11  state  vectors  for  day  251  and 
252  produced  RMS  position  differences  between  the  data  sets  of  about  1500  m.  One 
sigma  position  accuracy  of  54  m  in  the  vertical  direction  and  approximately  1400  m  in  the 
downtrack  direction  were  experienced.  Velocity  vector  magnitude  differences  during  this 
period  were  generally  +  Im/s,  with  a  RMS  velocity  difference  of  less  than  9  m/s.  One 
sigma  velocity  accuracies  of  approximately  4.2  m/s  in  the  vertical  direction,  2.3  m/s  in  the 
downtrack  direction  and  1.5  m/s  in  the  crosstrack  direction  were  experienced.  A  firm 
conclusion  regarding  GPS  accuracies  could  not  be  drawn  because  all  sources  of  error 
were  not  identified. 

GPS  position  and  velocity  data  was  generally  very  good;  however,  spurious  data 
points  were  present.  The  errors  included  in  these  data  points  were  sometimes  very 
significant,  and  outages  in  GPS  data  of  approximately  two  hours  appeared  to  occur  at 
least  once  per  day.  Despite  these  errors,  GPS  appeared  to  be  effective  in  producing  good 
state  vector  data  even  during  vehicle  maneuvers. 

Based  on  these  results  GPS  appears  to  be  an  excellent  navigation  source  for 
providing  Shuttle  state  vector  information.  Considering  that  this  data  can  be  obtained 
real-time  and  still  match  the  reference  so  closely,  real-time  GPS  state  vector  data  appears 
even  more  valuable.  An  important  caveat  must  be  applied  to  this  very  valuable  GPS  data: 
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another  navigation  source  such  as  INS  must  be  present  to  provide  a  check  against 
spurious  data  points  and  periods  of  outage. 

B.  RECOMMENDATIONS 

Analysis  of  state  vector  data  from  a  Precise  Positioning  Service  GPS  receiver 
should  be  conducted.  PPS  state  vector  data  should  provide  considerably  better  position 
and  velocity  data  which  is  required  during  rendezvous  and  proximity  operations.  STS-79 
will  provide  an  opportunity  to  assess  PPS  GPS  receiver  navigation  solutions. 

Before  integration  with  the  Shuttle  navigation  suite,  experimentation  with  a  simple 
filter  for  the  GPS  receiver  navigation  solutions  is  required.  Erroneous  GPS  state  vector 
data  was  sometimes  quite  severe.  A  process  such  as  a  Kalman  filter  is  required  to  remove 
navigation  solutions  outside  a  specified  tolerance. 

A  standard  library  of  code  in  a  commonly  accepted  language  such  as  Fortran 
should  be  adopted  for  frequently  required  astronautical  functions  such  as  conversion 
between  reference  frames.  In  particular,  a  simple  uniform  routine  for  transformation 
between  the  J2000  ECI  frame  and  the  WGS-84  ECEF  frame  using  general  precession 
according  to  IAU-1976  and  general  nutation  according  to  IAU-1980  without  considering 
pole  wander  or  the  irregular  rotation  of  the  Earth  is  required.  Algorithms  which  allow  real 
time  processing  of  data  such  as  GPS  navigation  solutions  should  be  emphasized.  In 
contrast,  current  references  emphasize  the  use  of  techniques  which  require  post-flight 
processing  due  to  their  requirement  of  celestial  observations. 

GPS  system  time  should  be  adopted  by  all  tracking  sources.  Ambiguity  caused  by 
UTC  being  rounded  to  the  nearest  whole  second  by  the  GPS  receiver  and  tracking  sources 
led  to  significant  errors,  particularly  in  the  downtrack  direction.  Since  GPS  system  time 
can  be  obtained  to  a  very  high  accuracy  by  both  the  Shuttle  receiver  and  a  receiver  at  the 
tracking  site,  GPS  time  would  provide  a  valuable  timescale  addition  to  all  tracking  data 
sets. 
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APPENDIX  A.  FORTRAN  PROGRAM  FOR  EDITING  PROPAGATED  FILE 
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100 


program  edtprop 


implicit  none  -  declare  variables 

double  precision  hours,xft,yft,zft,vxft,vyft,vzft 
double  precision  secs,x,y,z,vx,vy,vz 
integer  i,count 

open(5,j51e='prop.out',status-old')  -  open  datafile 

rewind(5) 

open(6,file='prop.m50', status-unknown')  -  open  new  data  file 

rewind(6) 


write(6,*)  1  -  convertTool  header  to 

transform  to  J2000 

count=0 

do  100  i=l,90000 

read  (5,*,end=200)  hours, xft,yft,zft,vxft,vyft,vzft  -  read  data 


secs=hours*3600.0 

x=xfl*0.3048 

y=yfl*0.3048 

z=zft*0.3048 

vx=vxft*0.3048 

vy=vyft*0.3048 

vz=vzft*0.3048 

write(6,l)  secs,x,y,z,vx,vy,vz 


-  covert  MET  hours  to  seconds 

-  convert  feet  to  meters 


-  write  to  new  data  file 


count=count+l 
100  continue 
200  continue 


close(5) 

close(6) 

1  format(fl6.6,6(2x,fl6.6)) 

end 
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APPENDIX  B.  FORTRAN  PROGRAM  FOR  EDITING  NAVG-11  STATE 

VECTOR  FILE 
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program  edtsv 


implicit  none  -  declare  variables 

double  precision  station,xft,yft,zfl,vxft,vyft,vzft 
double  precision  x,y,z,vx,vy,vz 
integer  sv,doy,hr,min,sec 

integer  year,dayoyear,month,dayomon,hour,minute,second 
integer  tmission,secs 
integer  i,count 

open(5,file='sv.edt',status='old')  -  open  data  file 

rewind(5) 

open(6,file='sv.m50',status- unknown')  -  open  new  data  file 

rewind(6) 

year=1995 
dayoyear=250 
month=09 
dayomon=07 
houi:=15 
minute=9 
second=0 

tmission=dayoyear*24*3600+hour*3 600+minute*60+second  -  MET  start 

write(*,  1)1  -  convertTool  header  to 

1  format(il)  tran^orm  to  J2000 

count=0 

do  100  i=l,90000 

read  (5,*,end=200)  station,sv,doy,hr,min,sec,  -  read  desired  data 

%  xft,yft,zft,vxft,vyft,vzft 

secs=(doy*24*3600+hr*3600+min*60+sec)-tmission  -  MET  conversion 

x=xft*0.3048  -  convert  feet  to  meters 

5P=yfl*0.3048 

z=zft*0.3048 

vx=vxft*0.3048 

vy=vyfl*0.3048 

v^=vzft*0.3048 

write(6,2)  secs,x,y,z,vx,vy,vz  -  write  to  new  data  file 


-  data  required  for  MET 
conversion 


count=count+l 
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100  continue 

200  continue 


cIose(5) 

close(6) 

2  fonnat(i7,6(2x,fl6.6)) 

end 
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APPENDIX  C.  STK  VEHICLE  FILE 
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stk.v.2.0 


BEGIN  Vehicle 

Name  STS69gpsdayl 

BEGIN  VehiclePath 


BEGIN  J4Perturbation 

EphemEpoch 

StartTime 

StopTime 

SemiMajorAxis 

AltitudeOfApogee 

AltitudeOfPerigee 

RadiusO£\pogee 

RadiusOfPerigee 

PeriodOK)rbit 

MeanMotion 

Inclination 

Eccentricity 

ArgOfPerigee 

LongAscenNode 

RightAscension 

TrueAnomaly 

MeanAnomaly 

TimePastAscenNode 

TimePastPerigee 

TimeStep 

OrbitGranularity 

NumberOfPasses 

SizeShapeType 

AscenNodeType 

AnomalyType 

EllipseType 

END  J4Perturbation 

BEGIN  PassDefh 

Break 

Latitude 

DisplayFlag 


7  Sep  1995  15:09:00.000000 

7  Sep  1995  15:09:00.00 
19  Sep  1995  00:00:00.00 
6748537.000000 
370400.000000 
370400.000000 
6748537.000000 
6748537.000000 
5517.284751 
0.001139 
28.500000 
0.000000000000 
0.000000 
-151.000000 
249.528639 
0.000000 
0.000000 
0.000000 
0.000000 
60.000000 
0.000000 
0 

AppPeriAlt 

Longitude 

TrueAnomaly 

Osculating 


Ascending 

0.000000 

Both 
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FirstPass 
RangeFirstPass 
RangeLastPass 
DisplayScheme 
ScenarioEpoch 
Passes 
END  Passes 
END  PassDefh 

END  VehiclePath 

BEGIN  Ephemeris 

NumberOfEphemerisPoints  17342 

ScenarioEpoch  1  Nov  1992  00.00:00.0 

EphemerisEcfTimePosVel 

89944473.000000  -6290411.999970  -10669.639545  2444894.749963  -944.961651  - 
7710.190409  -2382.491604 

89944485.000000  -6300850.499989  -92193.726556  2416095.999986  -847.046094  - 
7710.033323  -2420.407165 

89944522.000000  -6325375.499997  -331401.656250  2329124.749954  -591.753766 
7704.550020  -2514.460050 


1 

0 

2147483647 

AllEphemeris 

1  Nov  1992  00:00:00.0 


90020994.000000  -5868722.494942  883252.418377  3213828.661367  -1541.348168 

7851.489098  -89.026636 

90020994.000000  -5868722.494942  883252.418377  3213828.661367  -1541.348168 

7851.489098  -89.026636 

90020994.000000  -5868722.494942  883252.418377  3213828.661367  -1541.348168 

7851.489098  -89.026636 

END  Ephemeris 

BEGIN  Attitude 

BEGIN  ECFWLH 

AZIMUTH  0.000000 

END  ECFWLH 

END  Attitude 
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BEGIN  Swath 


ElevationAngle  0.000000 

HalfAngle  0.000000 

RepType  NoSwath 

END  Swath 


BEGIN  Constraints 


ConstraintMask 

0 

MinRange 

0.000000 

MaxRange 

1000000000.000000 

MinAzimuthAngle 

-90.000000  1 

MaxAzimuthAngle 

180.000000  1 

MinElevationAngle 

-45.000000 

MaxElevationAngle 

45.000000 

MinGrazingAngle 

0.000000 

MaxGrazingAngle 

45.000000 

! 

MinGrazingAltitude 

185200.000000 

MaxGrazingAltitude 

1852000.000000 

MinGndElevationAngle 

0.000000 

MaxGndElevationAngle 

90.000000 

MinSunElevationAngle 

-90.000000 

MaxSunElevationAngle 

90.000000 

LightingCondition 

5 

END  Constraints 

BEGIN  Attributes 

MarkerColor 

6 

GroundTrackColor 

6 

SwathCoIor 

6 

LineStyle 

3 

MarkerStyle 

4 

END  Attributes 

BEGIN  Graphics 


Inherit 

On 

ShowLabel 

On 

ShowGndTrack 

On 
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END  Graphics 

BEGIN  Extensions 

BEGIN  Desc 
END  Desc 

BEGIN  Group 
END  Group 

BEGIN  Eclipse 

Penumbra  Off 

PenumbraColor  7 

PenumbraLineStyle  0 
PenumbraLineWidth  1 


Sunlight  Off 

SunlightColor  7 

SunlightLineStyle  0 

SunlightLineWidth  1 

Umbra  Off 

UmbraColor  7 

UmbraLineStyle  0 

UmbraLineWidth  1 


END  Eclipse 

BEGIN  Exclusion 
END  Exclusion 

END  Extensions 

BEGIN  SubObjects 

Class  Sensor 

Horizon 


END  Class 
END  SubObjects 


END  Vehicle 


APPENDIX  D.  MATLAB  M-FILE  FOR  COMPARING  GPS  AND 

PROPAGATED  DATA 
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%  Program;  diffdayl.m 
%  Jim  Jones 
%  25  May  96 

%  This  program  displays  x,y  and  z  components  of  GPS  and  reference 
%  track  position  and  velocity  data  for  one  day  of  data.  It  also 
%  displays  the  magnitude  of  the  position  vector,  or  altitude, 

%  from  the  center  of  the  earth  and  the  magnitude  of  the  velocity 
%  vector. 

%  Load  the  GPS  data  file  for  one  day 
load  gpsdayl.j20  -ascii 
gpsin=gpsdayl; 

%  Create  time,  position  component  and  velocity  component  vectors. 

gpssec=gpsin(;,l); 

gpsx=gpsin(:,2); 

gpsy=gpsin(:,3); 

gpszr=gpsin(;,4); 

gpsxdot=gpsin(:,5); 

gpsydot=gpsin(:,6); 

gpszdot=gpsin(:,7); 

%  Load  the  reference  track  data  file  for  one  day 

load  propdayl.j20  -asdi 

propin=propdayl; 

%  Create  time,  position  component  and  velocity  component  vectors. 

propsec=propin(: ,  1 ); 
propx=propin(;  ,2); 
propy=propin(;,3); 
propz=propin(:  ,4); 
propxdot=propin(:  ,5); 
propydot=propin(:  ,6); 
propzdot=propin(:,7); 
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%  Plot  position  components  vs.  time 


figure(l) 

subplot(3,l,l) 

plot(gpssec,gpsx,*w*',propsec,propx,'w') 
xlabel('Day251  Mission  Seconds') 
ylabel('X  Position  (m)') 
grid 


subplot(3,l,2) 

plot(gpssec,gpsy,'w*',propsec,propy,'w') 
xlabelCDay  251  Mission  Seconds') 
ylabel('Y  Position  (m)') 
grid 

subplot(3,l,3) 

plot(gpssec,gpsz,'w*',propsec,propz,'w') 
xlabel('Day  251  Mission  Seconds') 
yIabeI('Z  Position  (m)') 
tol=-l; 

legend('GPS','Reference',toI) 

grid 

orient  tall 

%print_orig  -dgif8  251pos.gif 
print 

%  Plot  velocity  components  vs.  time 

figure(2) 

subplot(3,l,l) 

plot(gpssec,gpsxdot,'w*',propsec,propxdot,'w') 
xlabel('Day  251  Mission  Seconds') 
ylabelCX  Velocity  (m/s)') 
grid 

subplot(3,l,2) 

plot(gpssec,gpsydot,'w*',propsec,propydot,'w') 
xlabel('Day  251  Mission  Seconds') 
ylabel('Y  Velocity  (m/s)') 
grid 
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subplot(3,l,3) 

plot(gpssec,gpszdot,'w*',propsec,propzdot,'w') 
xlabelCDay  251  Mission  Seconds') 
ylabel('Z  Velocity  (m/s)') 
tol=-l; 

legend('GPS','Reference',tol) 

grid 

orient  tall 

%print_orig  -dgifS  251vel.gif 
print 

%  Calculate  position  vector  magnitudes  (altitude) 

gpsr=sqrt(gpsx.^2  +  gpsy.^2  +  gpsz.^2); 
propr=sqrt(propx.''^2  +  propy.'^2  +  propz.'^2); 

%  Calculate  velocity  vector  magnitude 

gpsv=sqrt(gpsxdot.''^2  +  gpsydot.^2  +  gpszdot.'^2); 
propv=sqrt(propxdot.^2  +  propydot.^2  +  propzdot.'^2); 


%  Plot  position  and  velocity  vector  magnitudes 

figure(3) 

subplot(2,l,l) 

plot(gpssec,gpsr,'w*',propsec,propr,'w') 

xlabelCDay  251  Mission  Seconds') 

ylabeI(Tosition  Vector  Magnitude  (Vehicle  Altitude)  (m)') 

tol=-l 

legend('GPS','Reference',tol) 

grid 

subplot(2,l,2) 

plot(gpssec,gpsv,’w*',propsec,propv,'w') 
xlabel(Day  25 1  Mission  Seconds') 
ylabel('Velocity  Vector  Magnitude  (m/s)') 
tol=-l 

legend('GPS','Reference',tol) 

grid 

orient  tall 

%print_orig  -dgifS  251niag.gif 
print 

%  Plot  close  up  of  position  data 
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figure(4) 

plot(gpssec,gpsx,'w*',propsec,propx,'wo') 
xIabel('Day  251  Mission  Seconds') 
ylabelCX  Position  (m)') 
tol=-l 

legend('GPS','Reference',tol) 

grid 

axis([(5e4+20)  (5e4+220)  -0.57e7  -0.47e7]) 

%print_orig  -dgifS  251close.gif 

print 


118 


APPENDIX  E.  MATLAB  M-FDLE  FOR  PLOTTING  GPS  ORBIT  IN  3-D 
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%  Program:  orbpltl.m 
%  Jim  Jones 
%  25  May  96 

%  This  program  creates  a  3D  plot  in  J2000  coordinates  of  GPS  navigation 
%  solutions  for  one  day  of  flight  data. 

%  Load  GPS  ephemeris  file  in  X,Y,Z  J2000  coordinates 

load  gpsday0.j20  -ascii 

out  =  gpsdayO; 

%  Plot  ephemeris  from  data  file 

xl  =  [0  le7]; 
x2  =  [0  0]; 
x3  -  [0  0]; 

yl  =  [0  0]; 
y2  =  [0  1e7]; 
y3  =  [0  0]; 

zl  =  [0  0]; 
z2  =  [0  0]; 
z3  =  [0  le7]; 

figure(l) 

plot3(xl,x2,x3,'w-’,yl,y2,y3>-',zl,z2,z3,'w-',out(;,2),out(:,3),out(;,4),'w.') 

title('GPS  Orbit  for  Day  250  in  J2000  Coordinates') 

xlabel('X  Axis  (m)') 

ylabel('Y  Axis  (my) 

zlabel('Z  Axis  (m)*) 

%print_orig  -dgi®  250orb.gif 
print 
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APPENDIX  F.  MATLAB  M-FILE  FOR  COMPARING  GPS  AND  NAVG-11 
STATE  VECTORS  WITH  MATCtONG  TIMES 
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%  Program:  difimatl.m 
%  Jim  Jones 
%  25  May  96 

%  This  program  plots  differences  between  GPS  navigation  solution 
%  data  and  Rockwell  state  vector  data  which  share  the  same  time. 
%  Both  data  sets  are  in  J2000  coordinates. 


%  Load  GPS  data  file 
load  gpsmatchl.j20  -ascii 
gpsin=gpsmatch  1 ; 

%  Create  time,  position  component  and  velocity  component  vectors. 

gpssec=gpsin(;,l); 

gpsx=gpsin(:,2X 

gpsy=gpsin(:,3); 

gpsz=gpsin(:,4); 

gpsxdot=gpsin(:,5); 

gpsydot=gpsin(:,6); 

gpszdot=gpsin(:,7); 

%  Load  Rockwell  state  vector  data 

load  svmatchl.j20  -ascii 

propin=svmatch  1 ; 

%  Create  time,  position  component  and  velocity  component  vectors. 

propsec=propin(: ,  1 ); 

propx=propin(:,2); 

propy=propin(;,3); 

propz=propin(:,4); 

propxdot=propin(:  ,5); 

propydot=propin(;  ,6); 

propzdot=propin(:,7); 
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%  Identify  sources  of  Rockwell  state  vectors 
%  (Sources  corresponding  to  indices  of  matrix  propin) 


%  Legend 

estr=[l  3  5  7  9  1 1  16  21  24  27  30  32  34  36];  %  * 

wstr=[2  4  6  8  10  13  17  18  22  28  35  37];  %  o 

bdqc=[12  15  20  26  31];  %x 

wlrc=[14  19  25];  %  + 

kmtc=[23  29  33];  %  • 


%  Plot  position  components  vs.  time 

figure(l) 

subplot(3,l,l) 

plot(gpssec,gpsx,'w*',propsec,propx,'wo') 
xlabel(T)ay  251  Mission  Seconds') 
ylabelCX  Position  (m)-) 
grid 

subpIot(3,l,2) 

plot(gpssec,gpsy,'w*',propsec,propy,'wo') 
xlabel(T)ay  251  Mission  Seconds') 
ylabel('Y  Position  (m)') 
grid 

subplot(3,l,3) 

plot(gpssec,gpsz,'w*',propsec,propz,'wo') 
xlabel('Day  251  Mission  Seconds') 
ylabel('Z  Position  (m)') 
tol=-l; 

legend('GPS','Reference',tol) 

grid 

orient  tall 

%print_orig  -dgjfS  251posm.gif 
print 

%  Plot  velocity  components  vs.  time 

figure(2) 

subplot(3,l,l) 

plot(gpssec,gpsxdot,'w*',propsec,propxdot,'wo') 

xlabel('Day  251  Mission  Seconds') 

ylabelCX  Velocity  (m/s)') 

grid 
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subplot(3,l,2) 

plot(gpssec,gpsydot,'w*',propsec,propydot,'wo') 
xlabelCDay  251  Mission  Seconds') 
ylabel(T  Velocity  (m/s)') 
grid 

subplot(3,l,3) 

plot(gpssec,gpszdot,'w*',propsec,propzdot,'wo') 
xlabelCDay  251  Mission  Seconds') 
ylabel('Z  Velocity  (m/s)') 
tol=-l; 

legend('GPS','Reference',tol) 

grid 

orient  tall 

%print_orig  -dgifS  251velm.gif 
print 

%  Calculate  position  vector  magnitudes  (altitude) 

gpsr=sqrt(gpsx.^2  +  gpsy.^2  +  gpsz.''2); 
propr=sqrt(propx.''^2  +  propy.^2  +  propz.'^2); 

%  Calculate  velocity  vector  magnitude 

gpsv=sqrt(gpsxdot.^2  +  gpsydot.'^2  +  gpszdot.^2); 
propv=sqrt(propxdot.''^2  +  propydot.^2  +  propzdot.''^2); 


%  Plot  position  vector  and  velocity  magnitudes 

figure(3) 

subplot(2,l,l) 

plot(gpssec,gpsr,'w*',propsec,propr,'wo') 

xlabel(Day  25 1  Mission  Seconds') 

ylabelCPosition  Vector  Magnitude  (Vehicle  Altitude)  (m)') 

grid 

subplot(2,l,2) 

plot(gpssec,gpsv,'w*',propsec,propv,'wo') 
xIabel(Day  25 1  Mission  Seconds') 
ylabel('Velocity  Vector  Magnitude  (m/s)') 
tol=-l 

legend('GPS','Reference',tol) 

grid 

orient  tall 

%print_orig  -dgifS  251magm.gif 
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print 

%  Calculate  position  component  differences 

posdifl6c=gpsx-propx; 

posdiflfy=gpsy-propy; 

posdifi^gpsz-propz; 

%  Plot  position  component  differences 

%  Data  source  legend 

%  estr-  * 

%  wstr-  o 
%  bdqc-  X 
%  wire  -  + 

%  kmte-  . 

figure(4) 

plot(propsec(estr),posdiffs:(estr),'w* • 

propsec(wstr),posdifiEx(wstr),'wo',propsec(bdqc),posdifBc(bdqc),'wx', . . . 
propsec(wlrc),posdiffx(wlrc),'w+',propsec(kmtc),posdifiBc(kmtc),'w.') 
xlabelCDay  251  Mission  Seconds') 
ylabelCX  Position  Difference  (m)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC','KMTC',tol) 

%print_orig  -dgifS  251difxm.gif 
print 

figure(5) 

plot(propsec(estr),posdiffy(estr),'w*',. . . 

propsec(wstr),posdiffy(wstr),'wo',propsec(bdqc),posdifiy(bdqc),'wx', . . . 
propsec(wlrc),posdiffy(wlrcVw+',propsec(kmtc),posdiflfy(kmtc),'w.') 
xlabel('Day  251  Mission  Seconds') 
ylabel('Y  Position  Difference  (m)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC','JCMTC',tol) 

%print_orig  -dgifS  251difym.gif 
print 
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figure(6) 

plot(propsec(estr),posdiffe(estr),'w*  V  •  • 

propsec(wstr),posdifiz(wstr),'wo',propsec(bdqc),posdiflfe(bdqc),'wx', . . . 
propsec(wlrc),posdiffe(wlrc),'w+',propsec(kmtc),posdifi&:(kmtc),'w.') 
xlabel('Day  251  Mission  Seconds') 
ylabel(’Z  Position  Difference  (m)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','\VLRC','KMTC',tol) 

%print_orig  -dgifS  251difem.gif 
print 


%  Calculate  velocity  component  differences 

veldifBc=gpsxdot-propxdot; 

veldiffy=gpsydot-propydot; 

veldififz=gpszdot-propzdot; 

%  Plot  velocity  component  differences 

figure(7) 

plot(propsec(estr),  veldifBc(estr),'w*', . . . 

propsec(wstr),veldifBc(wstr),'wo',propsec(bdqc),veldif&:(bdqc),'wx', . . . 
propsec(wlrc),veldifix(wlrc),'w+',propsec(kmtc),veldifBc(kmtc),'w.') 
xlabel(Day  251  Mission  Seconds') 
ylabelCX  Velocity  Difference  (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR’,'BDQC','WLRC',’KMTC',tol) 

%print_orig  -dgifS  251dfxdm.gif 
print 

figure(8) 

plot(propsec(estr),  veldiffy(estr),'w*', . . . 

propsec(wstr),veldiffy(wstr),*wo',propsec(bdqc),veldiffy(bdqc),'wx', . . . 
propsec(wlrc),veldiffy(wlrc),'w+',propsec(kmtc),veldiffy(kmtc),'w.') 
xlabel(Day  251  Mission  Seconds') 
ylabelCY  Velocity  Difference  (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR’,'BDQC,'WLRC','KMTC',tol) 

%print_orig  -d^fS  251dfydm.gif 
print 
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figure(9) 

pIot(propsec(estr),veIdifl&:(estr),'w*', . . . 

propsec(wstr),  veldiflEz(wstr),'wo',propsec(bdqc),  veldifEz(bdqc),'wx', . . . 
propsec(wlrc),veldiffi:(wlrc),'w+',propsec(kmtc),veldifiz(kmtc),'w.') 

xlabel('Day  251  Mission  Seconds') 
yIabel('Z  Velocity  Difference  (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC’,'KMTC',tol) 

%print_orig  -dgifS  251dfzdm.gif 
print 


%  Calculate  position  and  velocity  magnitude  differences 

diffr=gpsr-propr; 

diffv=gpsv-propv; 

%  Plot  position  and  velocity  magnitude  differences 
figure(lO) 

plot(propsec(estr),dif&(estr),'w*',. . . 

propsec(wstr),diffr(wstr),'wo',propsec(bdqc),difir(bdqc),'wx', . . . 
propsec(wlrc),dif&(wlrc),'w+',propsec(lantc),diflfr(kmtc),'w.') 

xlabel(Day  251  Mission  Seconds') 

ylabelCPosition  Vector  Magnitude  (Altitude)  Difference  (m)') 

grid 

tol=-l; 

legend('ESTR','WSTR','BDQC',’WLRC','KMTC',tol) 

%print_orig  -dgifS  251difrm.gif 
print 

figure(ll) 

plot(propsec(estr),difiv(estr),'w*',. . . 

propsec(wstr),difR^(wstr),'wo',propsec(bdqc),difiV(bdqc),'wx', . . 
propsec(wlrc),diffv(wlrc),'w+',propsec(kmtc),diffv(kmtc),'w.') 
xlabel(Day  251  Mission  Seconds') 
ylabel('Velocity  Vector  Magnitude  Difference  (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC','KMTC',tol) 

%print_orig  -dgifS  251difvm.gif 
print 
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%  Calculate  position  and  velocity  RMS  differences 
rmsr=sqrt(posdifEc/'2  +  posdifiy.^2  +  posdiffe.^2); 
misv=sqrt(veldiffic.'^2  +  veldiSy.^2  +  veldifiy.^2); 

%  Plot  position  and  velocity  RMS  differences 
figure(12) 

plot(propsec(estr),rmsr(estr),'w*',.  •  • 

propsec(wstr),nnsr(wstr),'wo',propsec(bdqc),rmsr(bdqc),'wx',. 
propsec(wlrc),nnsr(wlrc),'w+',propsec(knitc),rmsr(kmtc),'w.') 
xlabelCDay  251  Mission  Seconds') 
ylabeK’Position  RMS  Difference  (m)') 
grid 
tol=-l; 

iegend('ESTR’,'WSTR','BDQC','WLRC','KMTC',tol) 

%print_orig  -dgifB  251drmsm.gif 
print 

figure(13) 

plot(propsec(estr),rmsv(estr),'w*’, . . . 

propsec(wstr),misv(wstr),'wo',propsec(bdqc),rmsv(bdqc),'wx', 
propsec(wlrc),nnsv(wlrc),'w+',propsec(kmtc),rnisv(kmtc),'w.') 
xlabel(Day  251  Mission  Seconds') 
ylabel('Velocity  RMS  Difference  (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC','KMTC',tol) 

%print_orig  -dgifS  25 1  dvmsm.^ 
print 
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APPENDIX  G.  MATLAB  M-FILE  FOR  TRANSFORMING  MATCHING  GPS 
AND  NAVG-11  STATE  VECTORS  TO  UEN  FRAME  FOR  COMPARISON 
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%  Program;  uen251.m 
%  Jim  Jones 
%  25  June  96 

%  This  program  plots  differences  between  GPS  navigation  solution 
%  data  and  NAVG-1 1  state  vector  data  which  share  the  same  time. 
%  Both  data  sets  are  in  the  Shuttle  UEN  frame. 

%  Load  the  GPS  data  file  for  one  day 

load  gpsmatchl.j20  -ascii 

gpsin=gpsmatchl; 

%  Create  time,  position  component  and  velocity  component  vectors. 

gpssec=gpsin(:,l); 

gpsx=gpsin(:,2); 

gpsy=gpsin(:,3); 

gpsz=gpsin(:,4); 

gpsxdot=gpsin(:,5); 

gpsydot=gpsin(;,6); 

gpszdot=gpsin(:,7); 

%  Load  the  reference  track  data  file  for  one  day 

load  svmatchl.j20  -ascii 

svin=svmatchl; 

%  Create  time,  position  component  and  velocity  component  vectors. 

refsec=svin(:,l); 

refx=svin(;,2); 

refy=svin(:,3); 

refe=svin(:,4); 

refxdot=svin(:,5); 

refydot=svin(:,6); 

refedot=svin(:,7); 
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%  Identify  sources  of  NAVG-1 1  state  vectors 
%  (Sources  corresponding  to  indices  of  matrix  svin) 


%  Legend 

esti=[l  3  5  7  9  1 1  16  21  24  27  30  32  34  36];%  * 
wstr=[2  4  6  8  10  13  17  18  22  28];  %  o 

bdqc=[12  15  20  26  31];  %x 

wlrc=[14  19  25];  %  + 

kmtc=[23  29  33];  %  . 


for  i=l:length(gpsx) 

%  Elements  required  to  calculate  DCM  for  GPS  data 

gpsR=sqrt(gpsx(i)"'2+gpsy(i)^2+gpsz(i)^2); 

gpsgamma=atan2(gpsy(i),gpsx(i)); 

gpsd=sqrt(gpsx(i)^2+gpsy(i)^2); 

gpsdelta=atan2(gpsz(i),gpsd); 

%  DCM  components  for  GPS  data 

gpsdcm  1  l=cos(gpsdelta)*cos(gpsgamma); 
gpsdcml2=cos(gpsdelta)*sin(gpsgamma); 
gpsdcml  3=sin(gpsdelta); 

gpsdcm2 1  =sin(gpsdelta); 

gpsdcm22=cos(gpsgamma); 

gpsdcm23=0; 

gpsdcm31=-sin(gpsdelta)*cos(gpsgamma); 

gpsdcm32=-sin(gpsdelta)*sin(gpsgamma); 

gpsdcm33=cos(gpsdelta); 

%  DCM  rows  for  GPS  data 

gpsdcml=[gpsdcmll  gpsdcml2  gpsdcml3]; 
gpsdcm2=[gpsdcm21  gpsdcm22  gpsdcm23]; 
gpsdcm3=[gpsdcm3 1  gpsdcm32  gpsdcm33]; 

%  3x3  DCM  for  GPS  data 

gpsdcm=[gpsdcml;  gpsdcm2;  gpsdcm3]; 
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% 


Create  GPS  position  and  velocity  vectors 


gpspos=[gpsx(i);  gpsy(i);  gpsz(i)]; 
gpsvel=[gpsxdot(i);  gpsydot(i);  gpszdot(i)]; 

%  Perform  transformation  from  J2000  coordinates  to  UEN 

gpsuenp=gpsdcm*gpspos; 

gpsuenv=gpsdcm*gpsvel; 

%  Elements  required  to  calculate  DCM  for  reference  data 

refR=sqrt(re&(i)'^2+refy(i)^2+refe(i)'^2); 

refgamma=atan2(refy(i),refx(i)); 

refd=sqrt(refx(i)^2+refy(i)'^2); 

refdelta=atan2(re£z(i),refd); 

%  DCM  components  for  reference  data 

refdcm  1 1  =cos(refdelta)*cos(refgamma); 
refdcml2=cos(refdelta)*sin(refgamma); 
refdcm  1 3=sin(refdelta); 

refdcntZ  1  =sin(refdelta); 

refdcm22=cos(refgamma); 

refdcm23=0; 

refdcm31=-sin(refdelta)*cos(refgamma); 
refdcm32=-sin(refdelta)*sin(refgamma); 
refdcm3  3=cos(refdelta); 

%  DCM  rows  for  reference  data 

refdcml=[refdcml  1  refdcml2  refdcml3]; 
refdcm2=[refdcm21  refdcm22  refdcm23]; 
refdcm3=[refdcm31  refdcm32  refdcm33]; 

%  3x3  DCM  for  reference  data 

refdcm=[refdcml;  refdcm2;  refdcm3]; 

%  Create  reference  position  and  velocity  vectors 

re^os=[re6c(i);  refy(i);  refe(i)]; 
re^el=[refxdot(i);  refydot(i);  refzdot(i)]; 
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%  Perform  transformation  from  J2000  coordinates  to  UEN 


refuenp=refdcm*rei^os; 

refuenv=refdcm*reivel; 

%  Calculate  state  vector  differences 

difp=gpsuenp-refuenp; 

difv=gpsuenv-refuenv; 

vertp(i)=di^(l); 

dtrlq)(i)=di^(2); 

xtrkp(i)=difi)(3); 

vertv(i)=difV(l); 

dtrkv(i)=difv(2); 

xtrkv(i)=difv(3); 


end 

%  Calculate  average  differences 

vertpavg=mean(vertp); 

dtrkpavg=mean(dtrkp); 

xtrkpavg=mean(xtrkp); 

vertvavg=mean(vertv); 

dtrkvavg=mean(dtrkv); 

xtrkvavg=mean(xtrkv); 

%  Calculate  median  differences 

vertpmed=median(vertp); 

dtrkpmed=median(dtrkp); 

xtrkpmed=median(xtrkp); 

vertvmed=median(vertv); 

dtrkvmed=median(dtrkv); 

xtrkvmed==median(xtrkv); 
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%  Calculate  standard  deviations  of  differences 

vertpstd=std(vertp); 

dtrkpstd=std(dtrkp); 

xtrkpstd=std(xtrkp); 

vertvstd=std(vertv); 

dtrkvstd=std(dtrkv); 

xtrkvstd=std(xtrkv); 

%  Plot  position  component  differences 

figure(l) 

plot(gpssec(estr),vertp(estr),'w*',gpssec(wstr),vertp(wstrX'wo',.  •  • 

gpssec(bdqc),vertp(bdqc),'\vx',gpssec(wlrc),  vertp(wlrc),'w+', . . . 
gpssec(kmtc),vertp(kmtc),'w. ') 
xlabelCDay  251  Mission  Seconds') 
ylabel('Vetical  Position  Difference  (m)') 
grid 
tol=-l; 

legend('ESTR',’WSTR','BDQC','WLRC','KMTC',tol) 

print 

figure(2) 

plot(gpssec(estr),dtrkp(estr),'w*',gpssec(wstr),dtrkp(wstr),'wo',... 
gpssec(bdqc),dtrkp(bdqc),'wx',gpssec(wlrc),dtrkp(wlrc),'w+',. . . 
gpssec(kmtc),dtrkp(kmtc),'w. ') 
xlabel(Day  251  Mission  Seconds') 
ylabel(Downtrack  Position  Difference  (m)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC','KMTC',tol) 

print 

figure(3) 

plot(gpssec(estr),xtrkp(estr),'w*',gpssec(wstr),xtrkp(wstr),'wo',. . . 
gpssec(bdqc),xtrkp(bdqc),'wx',gpssec(wlrc),xtrkp(wlrc),'w+',. . . 
gpssec(kmtc),xtrkp(kmtc),'w.') 
xlabel(Day  25 1  Mission  Seconds') 
ylabel('Crosstrack  Position  Difference  (m)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC','KMTC',tol) 

print 
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%  Plot  velocity  component  differences 
figure(4) 

plot(gpssec(estr),vertv(estr),'w*',gpssec(wstr),vertv(wstr),'wo', . . . 
gpssec(bdqc),vertv(bdqc),'wx',gpssec(wlrc),vertv(wlrc),'w+',. . . 
gpssec(kmtc),vertv(kmtc),'w. ') 
xIabelCDay  25 1  Mission  Seconds') 
ylabel('Vetical  Velocity  Difference  (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR';BDQC’,'WLRC';KMTC',tol) 

print 

figure(5) 

plot(gpssec(estrXdtrkv(estr),'w*',gpssec(wstr),dtrkv(wstr),'wo', . . . 
gpssec(bdqc),dtrkv(bdqc),'wx',gpssec(wlrc),dtrkv(wlrc),'w+',. . . 
gpssec(kmtc),dtrkv(kmtc),'w.’) 
xlabel(Day  25 1  Mission  Seconds') 
ylabel(Downtrack  Velocity  Difference  (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC','KMTC',tol) 

print 

figure(6) 

plot(gpssec(estr),xtrkv(estr),'w*',gpssec(wstr),xtrkv(wstr),'wo',. . . 
gpssec(bdqc),xtrkv(bdqc),'wx',gpssec(wlrc),xtrkv(wlrc),'w+',. . . 
gpssec(lantc),xtrk:v(kmtc),'w.') 
xlabel('Day  251  Mission  Seconds') 
ylabel('Crosstrack  Velocity  Difference  (m/s)') 
grid 
tol=-l; 

legend('ESTR','WSTR','BDQC','WLRC','KMTC',tol) 

print 
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